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CHAPTER  I 


BACKGROUND 


1.1  INTRODUCTION: 

The  objective  of  this  research  program  was  to  characterize  three- 
dimensional  aspects  of  fatigue  crack  closure  relative  to  part-through  and 
through-thickness  flaws.  Personnel  who  contributed  to  the  three  year 
research  grant,  along  with  publications  and  presentations  resulting  from  this 
effort  are  summarized  in  Appendix  A. 

The  work  is  motivated  by  the  well  known  fatigue  crack  closure 
phenomenon  whereby  naturally  occurring  fatigue  cracks  do  not  open  in  a 
linear  elastic  manner.  Since  proposed  by  Elber  [1-3],  considerable  research 
has  focused  on  developing  the  closure  mechanism  for  fatigue  life  pre'diction 
schemes.  As  emphasized  in  reviews  by  Paris  [ij,  Banarjee  (5j,  and  Suresh 
and  Ritchie  6',  understanding  crack  closure  is  of  fundamental  importance 
for  describing  many  aspects  of  fatigue  crack  growth. 

One  explanation  for  fatigue  crack  closure  involves  the  plastically 
deformed  region  atu'ad  of  the  crack  tip.  As  the  crack  pro[)ag;ites  through 
the  successive  plastic  zones  ahead  of  its  tip,  a  plastic  wake  is  formed  which 
results  in  rr'sidual  deformations.  These  deformations  in  turn  hold  the  crack 
faces  closed  during  portions  of  positive  applied  lo.ad  cycle,  and  reduce  the 


effective  load  for  the  remainder  of  the  cycle.  Figure  1.1  schematically  shows 
the  development  of  the  crack  tip  plastic  zones  as  the  crack  propagates 
through  the  specimen  and  the  resulting  plastic  wake  containing  the  residual 
dehmmations.  Figure  1.2  shows  how  the  effective  stress  range  is  reduce(l  in  a 
typical  load  cycle  by  crack  closure.  Here  is  the  stress  intensity  factor 

necessary  to  open  the  crack  tip  and  ‘^Kg^eciive  actual  stress  intensity 

factor  range  during  which  the  crack  tip  is  completely  open.  Since  crack 


e.xtension  can  not  occur  unless  the  crack  tip  is  open,  only 


cifect;  ve 


contributes  to  further  crack  propagation.  Thus,  fatigue  crack  closure  is  one 
of  the  controlling  factors  in  crack  extension. 

It  is  well  known  that  the  crack  tip  plastic  zone  is  larger  at  the  free 
surfaces,  where  plane  stress  conditions  occur,  than  at  the  center  of  a  thick 
specimen  where  plane  strain  conditions  prevail  (7).  Figure  1.3  schematically 
shows  the  through-the-thickness  variation  of  the  plastic  zone  in  a  thick 
specinren.  This  through-thickness  plastic  zone  size  variation  has  been  used 
to  explain  phenomena  such  as  thickness  dependent  fracture  toughness, 
fracture  surface  appearance,  and  thickness  related  fatigue  crack  behavior. 
The  larger  pla.stic  zone  at  the  specimen  surface  would  imply  that  the  closure 
elfect  is  more  pronounced  at  the  surface  than  at  the  interior,  resulting  in  a 
slower  crack  growth  rate  at  the  free  surface.  This  through-thickness 
variation  in  crack  growth  rates  is  commonly  called  the  tuimelling  efl'ect. 
'I'lie  effect  of  tfu;  st.atc!  of  stress  on  plastic  zone  size,  and  the  r('sulting  fatigue 
crack  growth  has  been  demonstrat(!d  with  v.ariable  am[)litude  h);uling 
experiments,  where  thin  specimens  have  longer  crack  growth  lives  than  thick 
members  !H-  1 1|. 
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In  addition  to  crack  closure  associated  with  the  plastic  wake  behind  the 
crack  tip,  two  other  closure  mechanisms  have  been  proposed:  asperity 
induced  and  oxide  iiulucinl  closur*'.  'I'lio  asperity  induced  closure  model  ;  12- 
11:  states  that  cra<-k  surface  roughness  keeps  the  crack  faces  propped  optui 
under  zero  load.  'The  maximum  plastic  zone  size  in  this  model  is  smaller 
than  the  grain  size,  while  the  size  of  the  fracture  surface  roughness  is  on  the 
same  order  as  the  crack  tip  displacement.  In  order  to  meet  the  requirement 
for  a  small  plastic  zone  .-^ize,  as[)erity  induced  closure  is  generally  observed  at 
low  crack  growth  rates  (on  the  order  of  10  ^  rnm/cycle).  When  the  fracture 
surface  size  is  the  doinin.ant  factor,  the  crack  tends  to  grow  in  a  zig-zag,  out 
of  plane  path,  leading  to  signilicant  Mode  H  displacements  promoting 
asperity  induced  closure. 

In  the  oxide  induced  closure  mechanism  [lo-lhj,  the  formation  of  an 
oxide  layer  just  behind  the  crack  tip  prevents  the  crack  surfaces  from 
closing.  .\.s  in  the  asperity  induced  closure  mechanism,  the  thickness  of  the 
oxide  layer  is  comparable  to  the  crack  tip  displacements.  During  the  closing 
phase  of  the  load  cycle,  early  contact  occurs  between  the  two  crack  faces 
due  to  the  presence  of  the  oxide  layer,  resulting  once  more  in  the  closure 
phenomenon.  Oxiiie  induced  closure,  like  the  asperity  model,  has  also  been 
observed  at  low  crack  growth  rates.  Since  both  asperity  and  oxide  induced 
closure  mechanisms  kt'ep  tl\e  crack  faces  open  under  zero  load,  these  are 
sometimes  referre<l  (o  as  ".N'on-ciosure"  models. 


1.2  REVIEW  OF  PRIOR  WORK: 


This  section  reviews  techni<}ues  that  have  been  employed  to 
characterize  fatimie  crack  closure.  IJoih  nuinerical  and  experimental 
methods  are  brlelly  discussed. 


1.2.1  Numerical  Studies: 

.Newman  studied  crack  closure  in  a  center-cracked  panel  under  cyclic 
loading  with  a  two-dimensional,  non-linear,  finite  element  model  [17,.  In  this 
study  the  material  was  assumed  to  be  elastic-perfectly  plastic,  and  the 
model  was  composed  of  two-dimensional  constant-strain  triangular  elements. 
It  was  observed  that  the  elernent-mesh  size  near  the  crack  tip  inlluenced  the 
prediction  of  the  magnitude  of  crack  closvire  and  opening  loads.  By  choosing 
an  appropriate  iinite-elcment-mesh,  the  actual  experimental  crack  growth 
rate  could  be  simulated.  .Mthough  the  linite-element  method  may  work  well 
for  closure  predictions,  the  analysis  is  often  complicated  and  may  rcfjuire 
long  computation  times. 

Cherrnahini  18;  p(?rformed  a  three-dimensional  elastic-plastic  finite 
ehunent  an.alysis  of  a  straight  through  crack  under  cyclic  tensile  loading. 
His  analysis  did  not  allow  for  any  curvature  of  the  crack  front.  The  crack 
was  simply  extemh'd  .at  tlu;  maximum  load  of  each  cycle,  and  no  attempt 
was  made  to  c.alculate  th<’  .amount  of  cr.ack  growth.  'IIk.'  results  of 
(.dierrnahini’s  analysis  .appeared  to  agree  well  with  experimental  observation, 
and  also  gave  some!  n<'W  insight  into  the  stresses  that  devedop  behind  the 
crack  tip.  I'nfortunately,  this  type  of  an.alysis  requires  an  enormous  amount 


of  computer  time  even  though  the  program  was  vectorized  to  run  on  the 
VPS-32  computer  at  NASA  Langley  llesearch  Center. 

.Newman  ilOj  developed  a  two-dimensional  mo<lel  that  would  calculate 
crack  growth  for  a  center  cracked  panel.  'I'lie  nuxhd  was  based  on  the 
Dugdale  concept,  but  was  modified  to  leave  plastically  d(d’ormed  material  in 
the  crack  wake.  The  model  was  used  to  predict  crack  growth  in  center 
cracked  tension  specimens  of  2219-T851  aluminum  alloy  aluminum  alloy 
material.  Thirteen  tests  were  performed,  including  live  dilferont  types  of 
spectrum  loadings.  The  model’s  predictions  were  in  good  agreement  with 
the  experimental  results.  The  ratio  of  predicted  to  <'xperimental  lives  ranged 
from  0.66  to  1.18.  The  running  titncs  for  this  model  ranged  from  2  to  15 
minutes  on  a  CDC-6600  computer,  which  makes  this  method  much  more 
usable  than  finite  element  analysis. 

Fleck  120]  has  employed  the  two-dimensional  finite  element  analysis 
developial  in  reference  jlTj  to  study  crack  closure  undcjr  plane  strain 
conditions.  From  this  study,  it  was  noted  that  i)lasticit y-induced  closure 
under  plane  stress  conditions  is  (piite  different  from  the  plane  strain  closure 
behavior.  Under  plane  stress  conditions,  the  crack  faces  were  found  to  open 
and  close  in  a  continuous  manner  (i<‘:  crack  closi's  continuously  from  its  tip). 
Under  plane  strain  conditions,  a  discontinuous  (dosure  [)henomenon  was 
obsiTved,  wh('r(d)y  the  crack  first  closes  .at  a  location  near  the  cr.ack  tip,  and 
iH'Xt  closes  far  tx'hind  the  crack  ti[>,  leaving  .a  g.ap  where  no  (  losure  occurs. 


1.2.2  Experimental  Studies: 

Many  experimental  measurements  of  fatigue  crack  opening  have  l)e('n 
performed  on  various  typers  of  specimens,  employing  t('chni(iu('s  such  as  crack 
mouth  opening  displacement  (('MOD)  measurements,  strain  gages,  [)u-h-ro<l, 
etc.  Some  detailed  aspects  of  these  methods  are  discussed  below. 

The  CMOD  gage  [21-23]  measures  the  crack  mouth  opening 
displacement  from  a  clip  gage  mounted  across  the  mout'i  td’  the  precracklng 
notch.  .\  plot  of  displacement  versus  applied  load  is  obt. -lined,  and  tlie  lo;ui 
at  which  linear  elastic  behavior  begins  represmits  the  closuri'  load. 

The  strain  gage  measurement  techni<iuc  21-2.')]  involves  bonding  one  or 
more  strain  gages  at  various  locations  across  tl><>  crack  surface.  In  some 
cases  strain  gages  are  also  mounted  on  the  back  face  of  the  specimen.  The 
signal  from  the  strain  gages  are  then  recorded  as  a  function  of  applied  load, 
and  the  closure  load  is  again  determined  as  the  point  where  the  load  versus 
strain  record  becomes  linear. 

'fhe  ultrasonics  method  26-27|  measures  the  changing  acoustic 
resistance  of  a  specimen  as  the  crack  op<‘ns  or  closes.  Th<'  intensity  of  the 
ultrasonic  signal  reflected  from  the  fatigm;  crack  varies  depending  on  the 
.-vmount  of  crack  closure  jiri'sent.  In  this  technique,  an  ultrasonic 
transmitter  is  placed  on  the  top  of  the  cracked  member,  and  a  receiver  is 
placed  opposite  the  transmitter  on  the  bottom  of  the  specimen.  .'\s  before, 
the  received  signal  intensity  is  plotted  vi'rsus  lo;id,  ami  the  closure  load  is 
determined. 

'I'he  potenli.-il  dilfercuice  a{)proach  128-291  measures  the  electric.'il 
resistance  of  a  spc'cimen,  which  is  also  proportional  to  opening  of  the  crack. 
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In  this  instance  the  metal  specimen  acts  as  a  part  of  an  electrical  circuit.  A 
constant  current  supply  is  provided  across  the  specimen,  and  the  si'^iial 
obtained  from  potential  [irobes  pl.ai'c'd  on  both  sides  of  tlie  cr.'ick  i>  recorded 
,as  a  function  of  applied  load.  It  has  txa-n  obs<'r\'ed  in  -^ome  a p[ilica t  lo:.-- 
that  llie  received  signal  may  be  misled  by  the  ()r(;senc(;  of  .a  laser  ol' 
instilating  oxide  f)n  the  crack  faces  which  prevent  electrical  contact.  Otlu  r 
complicating  factors  may  include  tlie  change  of  electrical  properties  of  tlie 
material  in  th(^  crack  tip  \iehl  /.one. 

The  push-rod  disphicement  gage  .'}()- Ill  has  btam  usiai  to  determine  the 
closurt!  load  at  a  single  [)olnt  iii.^ide  the  spt-cimen.  h'or  tliis  method  a  pu.-h- 
rod  assembly  is  hastened  i<j  the  specimen  by  drilling  two  p.ar.allel  holes  just 
behind  tin-  f.atigiK'  crack  front.  The  rehative  di.s[>lacemeni  of  the  hole 
bottoms  is  me.asured  with  a  twin  cantilever  clip  gages  via  the  push-rods. 
The  closure  load  is  then  determined  by  locating  the  linear  point  on  tlte 
load- (lisp  ha  cement  curve. 

'I'he  h'.lber  g.age  tll’-tbl  metisures  the  crack  opening  at  a  certain  dist.ance 
from  the  tip  on  t  hi'  free  surftices  of  a  specimen.  Here  the  hilber  gage 
(modifn'd  cli[)  gtige)  is  mount<'ii  across  the  crack  plain'  and  the  closure  loaii 
is  determined  from  a  load-dis[)!.acement  [)lot  .as  l)efor('. 

Tlie  interferometric  displacement  gage  (1D(’>)  usi's  a  laser  to 

measure  ridative  displacement  between  two  shallow  rel’ectlng  indent  .a  t  ions  or 
grooves  located  across  the  crack.  Interference  fringe  [latterns  are  created  by 
the  ditfrai'teii  lasi-r  beams,  and  the  motion  of  these  fringes  r('[)resent  the 


crack  surface  displacement.  This  techni<pie  has  [trova'ii  to  lie  .an  ellectivc 
nn'thod  for  measuring  crack  surface  displ,ac(>ments  due  to  its  high  resolution 
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capability. 

Other  methods  to  obtain  the  closure  loads  include  direct-observation 
using  electron  microscopy  39',  and  a  vacuum  infiltration  techidKpie  10  .  0 

should  noted  luire  that  all  of  lhes(?  techniques  only  determine  the  closure 
load  at  the  specimen’s  surface,  or  at  a  single  point  inside  the  specimen,  and 
can  not  determine  the  complete  through-the-thickness  variation  of  closure. 
Furthermore,  since  the  closure  levels  observed  for  a  single  test  specimen  vary 
with  the  rneasurerntmt  location  [11-12',  the  results  obtained  from  the  above 
technifiue  cati  differ  significantly. 

While  the  above  techniques  reflect  an  average  (global)  or  single  point 
behavior  of  the  entire  specimen,  optical  interferometry  has  been  employeil  to 
determine  the  complete  through-the-thickness  variation  of  crack  closure  in 
transparent  polymer  specimens  [iS-dS].  These  optical  interference' 
experiments  have  shown  that  for  a  through-thickness  or  part-through  flaw, 
the  crack  opening  loads  vary  along  the  crack  front. 

1.3  OBJECTIVE  OF  CURRENT  RESEARCH: 

fhe  main  objective  of  the  current  research  is  to  deti'rmine  how 
through-the-thickness  variation  in  crack  tip  plasticity  affects  fatigue  crack 
o[>ening,  and  alters  fatigue  crack  growth  in  Ir.ansparent  mode!  materials. 

Both  experimental  and  numeric.al  .approaches  were  employed.  As 
descrif)('d  lati’r,  complete  three-dimensional  crack  opening  profiles  were 
me.asured  by  the  use  of  optical  interferometry  for  throiigh-thickne.ss  flaws, 
corner  crackl'd  holes,  and  surface  fl.aws.  'fhese  point  wise  measurements  art' 
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then  compared  to  global  (average)  measurements  obtained  from  CMOD  and 
back  face  strain  techniques.  Two  types  of  polymer  model  materials,  namely 
Polymethylmethacrylate  (PMMA)  and  Polycarbonate  (PC'),  are  included  in 
this  r(!search  prograni.  Although  this  project  did  no  testing  of  metals,  tin.- 
polymer  data  are  compared  to  results  of  renc-9r>  (a  nickel  based  alloy) 
testing  by  Ashbaugh  [49]  in  section  6.2.  Differences  in  fatigue  crack  growth 
mechanisms  of  metals  and  polymers  are  discussed  in  section  1.4  of  reference 


An  algorithm  was  also  developed  here  to  analyze  crack  growth  and 
crack  closure  of  a  surface  flaw.  This  algorithm  was  based  on  Newman’s 
model  for  center  cracked  plates,  but  was  modified  for  the  three-dimensional 
surface  crack  case.  The  model  was  used  to  analyze  the  surface  crack 
experiments  and  the  results  are  compared  to  the  experimental  values. 


Srlifniatlc  v'h'W  of  the  crack  lip  plastic  wake  lorrT.atifi 
the  crack  tij',  resulting  in  compressive  resiiiual  stress 
hold  the  crack  faces  closed  during  portions  of  p(isiii\i 
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CHAPTER  2 

EXPERIMENTAL  PROCEDURE 

2.1  OPTICAL  INTERFERENCE  TECHNIQUE: 

This  section  describes  the  procedures  used  for  the  experimental  portion 
of  the  research.  The  main  experimental  method  employed  for  the  current 
research  involved  optical  interference  measurements  of  crack  opening  in 
transparent  polymer  specimens.  Both  polymethylmethacrylate  (PMMA)  and 
polycarbonate  (PC)  model  materials  were  employed. 

Optical  interference  occurs  in  a  thin  transparent  wedge  when  the 
reflection  of  light  rays  from  the  top  and  the  bottom  of  the  faces  of  the 
wedge  have  different  path  lengths  [51-56].  When  a  crack  is  present  in  a 
transparent  material,  an  air  film  wedge  is  formed  between  the  two  crack 
surfaces  and  may  cause  optical  interference  to  take  place.  As  schematically 
presented  in  Figure  2.1,  some  light  waves  travel  through  the  transparent 
specimen  and  are  reflected  back  by  the  top  surface  of  the  crack,  whereas 
other  waves,  following  a  different  path,  penetrate  the  top  surface  and  are 
reflected  by  the  bottom  surface  of  the  crack.  This  difference  in  path  lengths 
(phase  difference)  causes  interference  fringes  to  form.  Each  fringe  represents 
a  locus  of  points  which  have  the  same  displacement  between  the  crack  faces. 
The  interference  fringes  form  an  alternating  pattern  between  dark  an<l  light 


shades  and  are  referred  to  as  destructive  and  constructive  fringes 
respectively.  If  the  wavelength  of  the  monochromatic  light  source  is  known, 
the  crack  surface  displacements  may  be  computed  using  either  of  the 
following  optics  equations  [o6j. 


For  destructive  interference  (dark  fringes): 


V  =  -  X 
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(2.1) 


For  constructive  interference  (light  fringes): 
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(2.2) 


Here  V  is  half  the  crack  surface  separation  at  a  specific  constructive  or 
destructive  fringe  location,  n  is  the  respective  fringe  order,  and  X  is  the 
wavelength  of  the  monochromatic  light  source.  Here  a  sodium  vapor  lamp  is 
employed  for  the  interference  study  with  a  wavelength  of  8.89  x  10~^  cm. 
The  O-order  dark  fringe  is  defined  here  as  the  first  destructive  fringe 
indicating  the  closed  portion  of  the  crack,  while  the  O-order  light  fringe 
corresponds  to  a  total  crack  separation  2V  =  X/4. 


2.2  THROUGH-THICKNESS  FLAW: 


The  through-thickness  Hawed  specimens  studied  here  consisted  of  edge 
cracked  beams  loaded  in  four-point  bending  as  shown  in  P'igure  2.2.  The 
test  specimens  were  38  mm  high  (VV)  and  178  mm  long  (L).  The  specimen 
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thickness  (B)  varied  from  25  mm  thick  to  5  mm  for  PMMA  and  from  25  mm 
to  10  mm  for  PC.  The  specimen  configuration  for  the  5  mrn  mernbors  was 
modified  to  avoid  iri.stability  on  the  four-point  bending  (i.xtures.  Here  the 
test  section  of  a  25  rnrn  tiiick  mernlx'r  was  reduced  to  rnm  with  siilficient 
care  being  taken  to  avoid  significant  stress  concentration  build  up  at  tlie 
test  section. 

A  total  of  five  different  specimen  thicknesses  for  PMMA  and  three  for 
PC  were  studied.  The  specimens  were  cut  from  a  single  slicet  of  the 
respective  material  and  each  test  member  contained  a  5  mm  deep  V^-notch 
located  midway  along  one  specimen  edge.  The  PM\L\  specimens  were 
annealed  at  100°C  and  the  PC  members  at  138°C  for  24  hours,  and  slowly 
cooled  to  room  temperature  to  reduce  potential  residual  stresses  due  to  the 
machining  process.  One  end  of  each  specimen  was  then  polished  to 
transparency  to  allow  crack  plane  observation.  Two  metal  tabs  were  glued 
to  each  specimen  at  the  starter  V-notch  to  hold  the  clip  gage  (MTS  Model 
632.03B-30)  at  the  crack  mouth  for  the  CMOD  measurements.  The  back 
face  strain  measurement  technique  involved  bonding  of  one  or  two  1000  ohm 
strain  gages  (Type  MA-06-250BK-10C)  on  the  top  side  of  the  four-point  bend 
specimens  (Figure  2.2). 

Cyclic  loads  (Haversine  function)  were  applied  in  a  pure-bend 
configuration  at  3Hz  for  PMMA  and  lllz  for  PC  specimens.  Thi'se  relatively 
low  test  frecjuencies  were  chosen  to  avoid  local  crack  tip  melting  due  to  the 
poor  heat  conductivity  of  the  polymers.  For  these  through-thickness  flaw 
experiments,  cracks  were  grown  under  constant  conditions,  and  constant 
stress  ratio  (F^)  was  maintained,  via  automated  computer  controlled  test 
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procedures  which  employed  compliance  measurements  of  crack  length. 
When  the  crack  grew  to  a  suitable  size,  interferometry,  CMOl),  and  BKS 
measvirements  were  performed  to  <letermine  the  closure  load.  'I'his  procedure 
was  re[jeat<'d  at  diHerent  crack  sizes  for  each  specinKui. 


2.3  CORNER  CRACKED  HOLES: 

Tiie  corner  cracked  hole  specimens  consisted  of  I’MNLV  rectangular 
plates,  with  a  hole  located  midway  along  the  specimen’s  length  as  shown  in 
Figure  2..3(a).  The  specimens  were  190  mm  long  (2h),  89  mm  wide  (2b),  and 
19  mm  thick  (T)  with  a  hole  diameter  (D)  of  19  mm.  For  the  present  study, 
the  corner  cracks  were  assumed  to  have  a  quarter-elliptic  configuration  as 
shown  in  F'igure  2.3(b).  Here  the  angle  is  the  elliptic  angle  commonly  used 
to  define  points  along  the  flaw  border  for  analysis  purposes  [57).  The  crack 
dimensions  a  and  c  represent  the  crack  length  along  the  hole  bore  and  the 
free  surface  locations  respectively. 

Cyclic  tensile  loads  were  applied  at  3Hz  to  the  PMMA  corner  cracked 
hole  specimens  through  grips  ghied  and  bolted  to  the  specimen  ends.  Two 
strain  gages  were  mounted  on  either  side  of  the  plate  and  monitored  to 
ensure  uniform  tensile  loading.  .\11  corner  cracks  were  grown  under  constant 
amplitude  load  conditions  with  a  fixed  stress  ratio  of  0.1.  The  crack  plane 
was  viewed  via  a  mirror  mounted  at  l.V’  on  the  transparemt  end,  and  was 
photographed  periixiically.  When  the  crack  achieved  a  suitable  length,  the 
cycling  was  stopped  and  a  set  of  optical  interference  fringe  patterns  were 
photographed  for  different  applied  tensile  loads.  This  procetlure  was 
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repeated  several  times  until  a  sufRciently  large  corner  flaw  appeared.  At 
that  time  the  loading  condition  was  changed  to  [)ure-t)endiiig,  keeping  the 
same  R  value  as  that  for  the  remote  tension  [iortioi\  of  the  expcritncni. 
Interference  fringe  patterns  were  again  pi'riodically  photographed  for 
increasing  crack  lengths  in  the  bending  configuration. 


2.4  SURFACE  FLAWED  PLATES: 

The  surface  flawed  members  consisted  of  PM\LV  rectangular  plates 
with  a  surface  flaw  starter  notch  located  midway  along  the  specimen’s 
length  as  shown  in  Figure  2.4  (a-b).  The  specimens  were  200  mm  long  (2h), 
102  mm  wide  (2b),  and  25  mm  thick  (t).  For  the  present  study,  the  surface 
flaws  were  assumed  to  have  a  semi-elliptic  configuration  as  shown  in  Figure 
2.4  (c)  with  the  elliptic  angle  O  used  to  define  the  points  along  the  flaw 
border.  The  crack  dimensions  a  and  c  represent  the  crack  depth  and  half 
the  free  surface  length  respectively.  For  the  CMOD  measurement  purposes, 
two  metal  tabs  were  glued  to  each  specimen  at  the  starter  notch  to  hold  the 
clip  gage  at  the  crack  mouth. 

The  surface  flawed  specimens  were  subjected  to  constant  amplitude 
remote  bending  cyclic  loads  with  a  R  value  of  0.1.  The  crack  plane  was 
periodically  photographed  through  the  transparent  end.  When  the  surface 
flaw  achieved  a  suitable  dimension,  the  cycling  was  stopped  and  a  set  of 
interference  fringe  patterns  were  photographed  for  increasing  applied 
bending  stress.  In  addition  to  the  optical  interference  measurements, 
readings  from  the  CMOD  gage  were  recorded  as  a  function  of  applied  load 
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Interferometry  orronq^ment 


CHAPTER  » 


SURFACE  CRACK  CLOSURE  MODEL 


'rhis  section  describes  the  numeric^!  a[)pro;ich  ernployc'd  to  predict 
three-dirnerisional  aspects  of  crack  closure.  'I’he  [irocednre  used  in  this 
analysis  for  the  calculation  of  fatigue  growth  of  a  surface  crack  is  basically 
the  same  as  Newman’s  method  for  a  center  cracked  panel  il9j.  His  method 
br(;ak.s  the  plastic  zone  and  plastic  wake  into  elements  (see  figure  .'hi)  which 
are  created  by  yielding  in  the  crack  tip  plastic  zone.  Ihifortunately  the 
tasks  are  greatly  complicated  by  the  three-dimensional  nature  of  the  surface 
flaw  problem,  and  by  the  absence  of  appropriate  crack  surface  displacement 
ecpiations.  Therefore  many  of  Newman’s  center  cracked  panel  equations  are 
used.  The  major  calculations  in  determining  closure  of  a  surface  Haw  are: 
d.l  Plastic  Zone  Size 

3.2  Deformation  Within  the  Plastic  Zone 
.'5.3  Contact  Stresses  at  Minimum  l.oad 
3.1  Crack  Opening  Stress  Intensity  Factor 
3..')  Fatigue  Crack  Growth 

A[)[)lication  of  the  method  and  ('(piations  used  by  Newman  for  the  center 
crack  [irobhun  are  (h'scribt'd  in  ch.apter  2  of  rt-ference  ."iS.  ('h.apter  3  of  that 
th('sis  describes  the  surface  crack  closure  model  more  explicitly  and  inclmh's 
more  details  than  are  given  in  this  report. 


3.1  PLASTIC  ZONE  SIZE; 


I'lie  [)l;ist'n‘  zone  siz(>  .iroutid  I  In’  pcriiihcry  of  ;i  surf.icc  crack  varies. 
'I'iiis  is  ihic  to  f'Ao  factors.  'I'lic  i!<  'j,rcc  of  coiist  r.aiiit  i-lianv,' s  from  pl.aiic 
strc.ss  .at.  t  lie  fn  c  urface  lo  ac.,;-  pl.aac  -tr.aii.  in  t!ic  interior,  ;iinl  i  tic  ^tre-' 
intensit }■  f.actor  warie.s.  .\ri  f  f.npirlc.il  stre-^s  intendty  f.actor  --oliilicHi  1^ 
avail.able  for  a  surrac('  cr.ack  in  hetiding  and  or  tanision.  It  was  devclopcai  by 
Newman  and  Rtijn  from  t  liree-dimensional  iinite  i.demerU  results,  .and  is 
given  below  in  its  general  form  .aO  (see  ligiirt'  d.2): 

K'l  =  (S-j.  f-  11  .^i()  ;  •-|.a/C>)|''  (a/t].(a/c).(>-/i>),  .  (.'blj 

Here  Q  i.s  the  eHi[)Lical  sli.ape  f.actor,  !■'  is  the  liound.ary  < orrection  factor  for 
tension  loading,  and  th('  product  of  H  and  !•'  is  tlie  bound.ary  correction 
factor  for  bending. 

For  this  research  the  plastic  zone  size  is  calculated  at  discrete  points 
around  the  crack  edgt?  by  the  following  steps  (see  figure  3.;?)  : 

1.  Calculate  K  at  Haw  locations  of  interest  by  ecpiation  3.1 

2.  Calculate  an  ecpiivalent  crack  length  (C,,,|)  by 
determining  the  length  from  the  cr.ack  origin  (point  O) 
to  the  point  of  interest 

3.  estimate  the  stress  state  coetlicient  (o) 

•1.  Solve  for  the  crack  tip  plastic  zone  siz«*  (/<)  by  using 
the  surf.ace  II, aw  v.alues  of  K,  ,an<]  i  in 

.Newman’s  center  cracked  panel  e<piations. 

These  ste[)s  .are  repe.ated  for  e.ach  point  of  interest  ( 1 ,2,;?.  1  ..a.t)  in  figure  3.1) 
along  the  cr.ack  periphery. 


3.2  DEFORMATION  WITHIN  THE  PLASTIC  ZONE: 


The  plastie  deformivtioiis  (1>|)  within  the  [)l:istic  zone  arc  calculated  in  a 
manner  similar  to  the  plastic  zone  size  calciilat  ioti.  '['he  siirfactt  lltiw  values 
of  K,  (  and  i  are  used  in  Ni-wm.an's  <'<-nter  crticked  panel  e(iu;ilions. 

I'he  transition  from  [)lane  stress  to  plam;  strain  was  basecl  on  closure 
experiments  performed  on  lltrough-cracked  I’MMA  specimens  of  various 
thicknesses  (25,  19,  13,  and  10  millitneters).  As  described  in  more  detail  in 
section  4.1,  optictil  interferometry  was  us<td  to  determine  the  opening  load  at 
points  along  the  crack  front.  'I'he  highest  opening  loads  occurred  at  the  free 
surface,  where  plane  stress  occurs.  The  opening  loads  for  all  thicknesses 
decreased  linearly  until  reaching  zero  at  a  distance  approximately  2.5 
millimeters  from  the  free  surface. 

In  this  model  it  is  assumed' that  plane  stre.ss  occurs  at  the  free  surface, 
plane  strain  occurs  when  the  distance  from  the  free  surface  is  greater  than 
0.1  inches,  and  the  stress  state  coefficient  alpha  is  linearly  interpolated  for 
distances  up  to  0.1  inches.  Since  PM\'L\  has  a  large  poisson’s  ratio 
(approximately  0.4),  this  allows  greater  stres-ses  to  develop  ahead  of  the 
crack  tip.  Therefore  alpha  was  ;i.s.signed  a  value  of  5  where  plane  strain 
occurred  [58].  In  ilOj  Newman  recommended  a  value  of  3  for  plane  strain, 
but  he  was  modeling  the  fadiavior  of  aluminum  which  has  a  poisson’s  ratio 
of  0.3. 

3.3  CONTACT  STRESSES  AT  MINIMUM  LOAD: 

'I'he  methotl  of  (h'termining  the  contact  stresses  in  the  (elastic  wakt? 
employed  for  this  study  utilizes  Newman’s  center  crack  displacement 
formula  along  six  radial  lines  which  •'onnect  [)oint  O  in  figure  3.3  with  crack 
perimeter  points  1,2, 3, 4, 5,  and  (5.  'I'his  is  done  in  the  following  manner: 


1.  The  stress  intensity  factor  at  point  1  is 
calculated  for  the  minimum  load  condition  using 
(•([uat  ion  d. ! . 

2.  I'sing  tli(!  ('([uivalent  crack  U'ngth  (uilculaled 

earlier,  and  the  stress  inleuisily  factor  determined  hy 
step  1,  the  crack  surface  displacements  due  to  the  remote 
load  are  calculated  by  using  the  center  cracked  panel 
ecpiations. 

An  iterative  melliod  is  used  to  determine  the  element  contact 
stresses.  I'hese  stresses  cause  the  crack  surfaces  to  open  up 
to  the  element  lengths.  When  necessary,  the  strc^sses  are 
modified  to  account  for  yielding  and  crack  separation. 

4.  Proceed  to  the  next  radial  line  (2, 3, 4, 5, 6)  and  repeat  until  the 
contact  stresses  along  all  lines  have  been  calculated. 

Using  this  method  appears  to  give  reasonable  contact  stresses  near  the 
crack  edge.  The  stresses  near  the  surface  crack  origin  (point  O)  are  not 
expected  to  be  accurate,  however,  since  shear  stresses  will  cause  the 
displacement  of  a  slice  to  be  strongly  alTected  by  the  displacements  of  its 
neighboring  slices.  This  affect  is  less  significant  at  the  crack  edge  because 
the  slices  are  farther  apart.  Inaccurate  stresses  near  the  crack  origin  are 
not  expected  to  affect  overall  accuracy  in  most  cases.  Newman  reported 
that  only  elements  near  the  crack  tip  carried  significant  loads  in  his  center 
cracked  panels  [I9j  for  R  ratios  greater  than  zero.  (Of  course,  his  results 
were  restricted  to  througli-lhe-thickness  flaws.)  ,\lso  stresses  away  from  the 
crack  edge  have  a  smaller  elfect  on  the  crack  opening  stress  intensity  factor 
than  stresses  near  the  crack  edge. 


3.4  CRACK  OPENING  STRESS  INTENSITY  FACTOR: 


When  llie  element  stresses  have  been  determined,  the  craek  <)[)enin'^ 
str('ss  i[it.(>nsity  factor  ('an  then  be  ealcnlated.  liy  usim^  the  Nvt'iyht  ftmciion 
iiK'thod,  Matthc'ck,  et  al  .(jOj  devel()[)('d  an  aittoritlmi  to  del  ermine  liic  -o  n  ss 
intensity  factors  at  points  a  and  c  (K^  and  K^)  due  to  arbitrary  stresses  on 
the  crack  surface.  Their  equations  are  given  below  and  the  coordinates  art* 
shown  in  figure  3.4. 


l<c 


■Ih/K^.-cl 

;ih/Kr,-aj  ^^^^’^W(x,y)l''d,F^/,  tcidxdy 


(3.-J) 


Here  h  is  a  matt^rial  constant  equal  to  E  for  plane  stress,  and  E/(l-/'')  for 
plam;  strain.  The  .str(‘ss  intensity  factors  at  A  and  C  due  to  the  referemat 
case  load  tire  and  K^,..  The  partial  derivatives  (tLy/i’a  and  ('U^/f'c  are 
the  changes  in  crack  surface  displacement  (U^)  due  to  the  reference  load 
with  respect  to  changes  in  the  crack  dimensions  a  and  c.  The  stress 
distribution  T'iew(X'y)  's  the  new  loading  for  which  and  K(,  will  be 
determined.  The  loading  must  be  symmetric  ((T^„^,(x,y)  =  (Tj^^^(x,-y)). 
Equations  3.2  and  3.3  are  evaluated  by  the  following  steps: 

1.  Calculate  a  and  c  at  element  centers 

2.  Break  crack  surface  into  sections  (figure  3.5) 

3.  .\verage  ''ly/<’a,  ''lJ,./t'c  and  contact  stri'ss  for  each  si'ction. 

4.  Numerically  evaluate  integrals  for  each  section. 

5.  Sum  integrals  and  multiply  by  constants. 


( l.'ilculating  the  parti.'il  derivatives  is  the  most  diHicult  step  her('.  This 
is  accomplished  by  calculating  the  displacements  for  the  actual  crack 
dimensions  (a,c),  a  slightly  deeper  crack  (a  f  'a,  c)  and  .'i  sliglitly  wider 
crack  (a,  c  f  i'c).  'fhe  derivatives  are  set  ecpial  to  the  change  in 
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displacement  divided  by  the  change  in  crack  size.  The  procedure  for 
approximating  the  displacements  is  rather  complex,  and  the  details  are  given 
in  ref('rence  ,()()'.  'I'ln'  surf.ace  crack  is  assumed  to  di^pl.ace  like  an  edge 
crai’k  along  y  0,  and  like  a  cenU'r  cr.ack  .along  \  0.  'Tin'  cr.ack  o[>eni!i'g 

a  [)proximat  ion  developcnl  by  I’elroski  and  .\<'henl).acli  (51  tV)r  tuo- 
diincnsional  problems,  was  modified  by  .Mattheck  et  al.  jGOj  Tor  the  surface 
Haw  problem.  In  this  study  an  outer  fiber  bending  stress  of  100  psi  was  used 
as  the  reference  load.  'I'his  valu<>  is  used  internally  by  the  algorithm  to 
calculate  displacements.  It  has  lu)  elfect  on  the  values  of  K.^  and  K,.  in 
equations  .3.2  and  3.3.  The  displacement  (UJ  and  the  reference  stri.’ss 
intensity  factors  and  K^^.)  are  directly  pro})ortional  to  the  refercmce 

stress  ('^new(^>y))i  effects  cat, cel.  'I'hc  type  of  refereiun'  case  stress 

(bending  or  tension)  does  have  some  elfect  on  the  algorithm  and  this  is 
discussed  in  section  5.1..^).  The  crack  is  broken  into  sections  by  connecting 
wake  elements  to  their  neighbors  in  adjacent  slices.  The  stresses  and  partial 
derivatives  of  the  two  elements  are  averaged,  and  these  values  are  assumed 
to  be  constant  over  the  entire  section.  'I'his  allows  easy  evaluation  of  the 
double  integral  for  each  section.  'I'lie  stress  intensity  factor  to  cause  cr.ack 
opening  is  equal  to  the  minimum  K  plus  the  K  caused  by  the  contact 
stresses,  and  is  given  by  equations  3.1  and  3.,‘>  for  crack  locations  .\  aiul  C. 

KOP.,  =  K„.,„(at  a)  +  \M\(at  a  cau.'^ed  by  contact  strc'sses)  (if-  f) 

KOP(.  =  K^in(at  c)  +  NjK(at  c  caused  by  contact  stresses)  (-3.')) 

'I'his  giv('s  two  opening  stress  intensity  factors,  KOP^  .and  K()I\,.  Ttu'se 
values  .are  uscnl  to  calculate  the  crack  growth. 

'I'he  algorithm  for  det{'rmining  K  dm;  to  cont.act  stressc's  was  <'valuated 
by  applying  bending  and  tensile  stresses  over  the  cr.ack  surfaces.  I’ln'  results 
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from  the  algorithm  were  then  compared  to  the  Nevvmari-Kaju  solutions. 
Although  the  Newman  and  llaju  results  are  for  either  reiiiot*'  leir  ion  or 
bt'uding,  tluty  may  he  used  to  evaluate  the  conlaet  o  re'-s  olulions  wla  ti  'he 
crtiek  ftiees  are  loaded  with  a  uniform  teu-ion  or  linearis’  sarsinv  hi  :,.;':  ..’ 
stre.s.s.  This  wtis  done  for  various  erack  sti.ape'-  .and  -'i/e.s  (ere  I  aide  !>.  1  i. 
Typical  results  are  shown  in  ligtires  .'{.6  and  it. 7.  'I'liere  is  good  agreement 
for  small  cracks,  but  the  method  is  inadcijuale  for  large  ;i/t  ratios  (above 
0.3).  This  observation  agrees  with  .\Iatthe<-k’s  results  50.  There  is  slightly 
better  agreement  at  (joint  .\  tlnan  .at  the  fr< c  surl'toa  .  The  largest  a/l  ratio 
used  in  the  calculations  shown  in  this  r<'[>ori  is  0. .')(),  w  hitdi  is  marked  on  the 
figures  by  a  dotted  line.  It  sh(jid<l  be  [jointeii  out  th.it  whih'  these  figures 
provide  a  comparison  of  the  tilgorithtn  to  an  ;iccc]>t ed  ('tti[)iric;il  solution,  the 
loading  is  very  dilferetit  from  what  the  crack  f.aces  would  experience  from 
surface  contact  due  to  closure.  Therefortt  these  figures  may  not  necessarily 
reflect  the  accuracy  of  the  algorithm  due  to  actual  loads. 

3.5  FATIGUE  CRACK  GROWTH: 

Newman  Mt),  employt'd  in  his  model  a  rather  complex  etjuation  for  crack 
growth.  'I'his  e(iuation  utilized  five  couslants  that  were  needed  to  model 
crack  growth  at  near  threshoUl  ami  near  fracture  conditions.  Since  this 
project  did  not  analyze.'  these  silu;itions,  the  slfuph'r  crack  growth  law 
(jroposed  by  l']lb(;r  j3j  (('(pialion  3.0)  should  be  sullicient. 

(la/du  ==  AIAK..,,)"  (3.0) 

I'ilhi'r's  eepiation  (3.0)  is  int<gr.it«'d  to  <icti'rmine  tin;  [lumber  of  cycles 
re<iuired  for  the  crack  to  grow  ,'>  [jercenl  of  the  [jlaslic  zone  size  at  (joint  0 
(figure  3.3). 


I*, 


N  =  /  '  {da/lA(AK„r)"l}  (3.7) 

a 

In  ('(Illation  .3. .3,  AK’,,(f  ^  point  0)  -  While  tlie  crack  i'^ 

‘;rowinj;  at  p(jint  (i.  it  i.s  also  griiwing  at  all  the  other  points  along  the  crack 
edge.  It  i.s  assumed  tliat  the  crack  will  retain  it.s  .semiellipt  ical  sha[)e  at  all 
times,  ’rherefore  as  the  crack  grows,  it  is  only  necessary  to  calculat('  new 
values  of  a  and  c.  W'liile  the  crack  grows  5  percent  of  its  [ilastic  /one  si/.('  at 
point  t),  its  growth  at  point  1  is  ca''”  'ated  by: 

Ac  =  N  ^  =  NA(AK„r'/‘  (3.s] 

dn 

'riie  .\  in  eii'uation  .3.8  is  the  N  (cycles)  that  wa.s  calculated  in  e(|iiation  .3.7. 
'I'he  AK,  -  K|,,_.^,j(at  point  1)  -  KOP^.  AHu'r  these  calculations  are 
coniplet(',  the  crack  is  extended  (a  —  a  !-  .05/’  and  c  --  c  '  Ac)  and  steps 
3.1  through  3.5  are  repeated. 
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I'  igure  .'1.6:  Test  of  closure  K  algorithm  at  point  A, 

maximum  bending  stress=1000  psi,  a/c  =  2/3 
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I'igiire  3.7:  Test  of  closure  K  .-rlgorithm  at  point  C, 

maximum  bending  stress— 1000  psi,  a/c— 2/3 


■"Iv  tn  AaBllh  I’lil 


*<  .V 


1 


CHAPTER  4 


EXPERIMENTAL  RESULTS 


This  chapter  describes  fatigue  crack  closure  results  obtained  from 
various  measurement  techniques  for  four  different  crack  types.  These  crack 
types  are: 

l.l  Through  cracks  in  PJylXlA 

•1.2  Through  cracks  in  polycarbonate 

4.3  Corner  cracked  holes 

4.4  Surface  flaws  in  bending 

These  experiments  are  described  in  greater  detail  in  chapters  4,5,6,  and  7  of 
reference  oOh 


4.1  THROUGH  CRACKS  IN  PMMA: 

The  test  matrix  for  through  thickness  flaws  is  presented  in  Table  4.1. 
Sections  4.1.1  and  4.1.2  will  focus  on  one  particular  test  (PX-6)  while  section 
4.1.3  examines  the  effect  of  varying  thickness,  R-ratio,  and  on  closure. 

F.xperimcnt  PX-6  was  conducted  at  a  cyclic  baseline  stress  intensity 
(^max)  ^25  KPa-rn^^^.  When  the  crack  length  reached  11.6  mm,  the  cyclic 
stress  intensity  was  reduced  to  550  KPa-m^^^.  The  stress  ratio  (R  —  Kmin  / 
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^^max)  kept  constant  at  0.1  for  the  entire  experiment. 

I'^i^ure  t.l  prc'sents  varions  measures  of  crack  length  versus  elapsed 
cych's.  /Vv(‘r;ige  f  hrougli-f  fiickness  crack  hmgtlis  obtained  from  both 
complian<‘e  and  photogra [)hic  measurements  are  shown  along  with  the  mid- 
plane  and  free  surface  crack  dimensions  obtained  from  the  photographs. 
Note  here  that  the  drop  in  from  825  to  550  KPa-m'^^  when  the  crack 

length  v/as  11.6  mm,  gave  a  small  crack  retardation  region  (marked  by  N,i) 
followed  by  steady  state  crack  growth  at  of  550  KPa-  m‘'^^.  .Vlthough 

the  compliance  (through-thickness)  crack  length  is  not  significantly  affected 
by  the  reduction  in  the  cyclic  stress  intensity  value,  it  can  be  seen  from 
Figure  l.l  that  the  free  surface  crack  dimension  exhibits  more  retardation 
than  the  mid-plane  crack  lengtfi.  This  crack  tunnelling  behavior,  where  the 
crack  grows  faster  in  the  mid-plane  than  the  free  surface  location,  was  also 
observed  in  PC  specimens  subjected  to  tensile  overloads  [44].  For  the 
present  experiment,  crack  tunnelling  continued  well  past  the  retardation 
period. 


4.1.1  Interferometric  Measurements: 

Figure  4.2  presents  a  set  of  interference  fringe  patterns  photographed 
during  the  steady  state  crack  growth  at  an  average  (compli<ance)  crack 
length  of  10.3  mm.  Recall  that  each  fringe  represents  a  locus  of  points  of 
constant  crack  face  displacement.  Figure  4.2(a)  shows  the  crack  plane 
under  zero  load  and  clearly  indicates  a  residual  crack  opening  displacement 
field.  Note  here  that  the  fringe  pattern  occupies  the  central  portion  of  the 
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crack  plane  and  a  part  of  the  intc'rior  crack  front.  Thus,  that  particular 
section  of  the  interior  crai-k  plane  is  o|)ened  under  zero  load,  while  the 
region  near  the  fret-  snrf.ice  remains  (•lose<l.  As  the  a[)plic;d  load  is  increasc^d. 
the  O-order  (outer-most)  light  fringe  reaches  the  fret^  surface  crack  li[)  at  a 
value  of  0.3,  wlumr!  is  the  peak  value  of  K  during  the  steady 

state  cycling.  The  load  recpiired  for  this  outer-most  light  fringe  to  reach  the 
crack  tip  is  one  measure  of  the  load  required  for  the  crack  surfaces  to 
separate,  and  is  dcliiied  here  as  Kqi  (iiiLerferometric  opening  load).  I'he  Ivq[ 
/  ^max  value  measured  at  the  free  surfa<’e  was  0.3,  while  this  value  is  zero 
for  the  mid-plane  location,  where  the  crack  surfaces  are  propped  open  at 
zero  load. 

Figure  1.3  shows  the  open  crack  perimeter  as  a  function  of  applied  load. 
Recall  that  in  Figure  4.2  (a-d)  each  photograph  showed  the  opened  crack 
area  for  a  particular  applied  load.  Figure  4.3  is  obtained  by  overlapping  the 
opened  areas  of  the  crack  for  increasing  applied  load  for  a  given  crack  size. 
For  example,  the  area  enclosed  between  the  crack  front  and  the  boundary 
line  of  K  /  represents  the  opened  section  of  the  crack  under  an 

applied  K  of  0.04K,nax-  Similarly  the  area  enclosed  between  the  crack  front 
and  any  other  boundary  line  represents  the  opened  section  of  the  crack 
under  the  respective  applied  K  value.  Note  here  that  under  a  small  load, 
portions  of  the  crack  faces  are  closed  at  the  crack  front  and  along  the 
starter  notch.  Further  applied  load  eventually  separates  the  crack  front  at 
the  free  surfaces. 

Crack  opening  displacements  are  examined  in  more  detail  by  digitized 
measurements  of  the  interference  fringe  patterns  from  Figure  4.2.  One 
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means  for  presenting  the  crack  opening  profile  is  to  plot  the  fringe  order  as  a 
function  of  distance  from  th<“  crack  tip  for  diliVrent  applied  loads,  as  seen  in 
h'igurt'  1.4  and  4. a,  giving  the  crack  suri‘a<-('  s(;paratlon  as  a  fiinclitMi  ol’ 
position.  The  displacements  are  ('X|)res.se<l  lu're  in  I'rlngc  order  units, 
although  other  dimensions  may  be  obtaiiuui  from  liquation  2.2  for  light 
fringes.  Two  displacement  profiles  were  obtained  for  each  load  sequence; 
one  measured  at  the  specimen’s  rnid-plane,  and  another  crack  opening 
profile  obtained  at  the  free  surface.  In  these  crack  o[)ening  proliles,  the 
physical  (original)  crack  tip  serves  as  the  plot  origin. 

Figure  4.4  presents  the  crack  opening  profiles  rneasurerl  at  the 
specimen’s  mid-plane  (interior)  for  an  average  crack  siz(;  of  10.3  mm.  .Since 
it  was  difficult  to  determine  the  fringe  number  near  the  mid-platie  crack  tip 
due  to  close  fringe  spacing,  as  seen  from  Figure  4.2,  the  fringe  order  near 
this  mid-plane  crack  tip  was  established  by  counting  the  fringe  number  from 
the  specimen’s  free  surface  location.  Each  curve  in  Figure  4.4  represents  a 
different  applied  load,  and  gives  the  total  separation  (2V)  between  the  two 
crack  surfaces  as  a  function  of  distance  from  the  crack  tip.  The  load  which 
causes  the  curve  to  pass  through  the  origin,  giving  complete  crack  tip 
separation,  is  referred  to  here  as  the  interferometric  opening  load  (K()|)  for 
the  particular  crack  location.  In  Figure  4.4,  the  crack  tip  in  the  specimen’s 
interior  is  separated  under  zero  applied  load.  I’igure  I..’)  shows  the 
corresponding  crack  opening  profiles  measured  at  the  s[)ecimen’s  frt'c 
surface.  .Note  here  that  crack  surfaces  are  closed  at  the  free  surface  until 
the  applied  K  equals  271  KPa-m*'^^  (ie:  Kqi  /  Emix  ‘  ^  Crack  mouth 
opening  displacement  for  an  applied  K  value  of  271  KFa-rn’'^^,  obtained  from 
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a  clip  gage  mounted  at  the  crack  mouth  is  shown  in  Figure  -I..').  The 
presence  of  the  start('r  notch  and  rough  precracking  surfaces  prc'veriied  the 
measurement  ol'  crack  mouth  opening  via  o[)iical  interfer(  ace.  In  addiiion, 
('lastic  I'rack  mouth  separation  computed  from  a  finite  elem.enl,  a.nal\  -l>  •’>() 
with  both  plain;  stress  and  plane  strain  assumptions  are  presented  in  l  igurt; 
1.5  along  with  the  clip  gage  reading.  Here  the  measurements  obtained  from 
the  clip  gage  and  Unite  element  analysis  are  converted  to  fringe  order  units. 

■Note  from  Figure  4.4  that  the  crack  surface  separation  is  a  in:'\itiuim 
at  a  certain  distance  from  the  crack  tip.  'I'his  maximum  crack  surface 
separation  location  moves  further  away  from  the  crack  tip  for  increasing 
load  as  seen  from  Figure  4.4.  Under  a  sufiiciently  high  load,  this  'bump’, 
where  the  crack  surface  separation  is  a  maximum,  moves  to  the  crack 
mouth.  Thus  under  a  large  load,  the  crack  surface  separation  only  increases 
as  one  moves  further  away  from  the  crack  tip.  Some  possible  factors 
contributing  to  the  fact  that  crack  separation  does  not  attain  its  maximum 
value  at  the  crack  mouth  under  ’low’  loads,  may  be  the  high  closure  stresses 
present  at  the  precracking  region  and  the  effect  of  free  surface  closure  on 
the  mid-plane  opening  behavior.  Figure  4.5  shows  that  unlike  the  crack 
opening  profile  at  the  mid-plane,  crack  surface  separation  reaches  a 
maximum  at  the  crack  mouth.  Furthermore,  the  magnitude  of  crack 
separation  is  much  smaller  for  the  free  surface  location  than  that  measured 
at  the  rnid-plane  region.  This  dilference  in  crack  opening  is  [iresumably  due 
to  the  larger  plastic  zone  at  the  specimen’s  free  surface. 

Figure  4.6  presents  the  interferometric  crack  tip  opening  load 
determined  at  different  locations  through  the  specimen  thickness  for  an 
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average  (compliance)  crack  length  of  10.3  mm.  Here  Kqi  values 
normalized  with  the  rnaxinnim  cyclic  stress  intensity  and  distancf's 

from  the  specimen's  free  surf.-icc  arc  norm.alized  with  the  ‘'jx'cimcn  thickiic,s 
(B).  The  abscissa  origin  r(‘pr(-senls  the  s[)ecin;en's  free  surface  locatii;n. 
while  x/B  ().")  specilies  t  lie  mid-plane  (interior)  location.  .Note  from 
Figure  4.6  that  the  opening  load  decrt.ases  rapiiily  as  one  movi's  into  the 
specimen  thickness.  .Moreover  about  70%  of  the  specimen’s  thickness  is  open 
under  zero  load.  One  interpretation  of  Figure  1.6  may  .suggest  that  dO'  c  of 
the  specimen’s  thickness  (l.'i'  c  from  (  ach  free  surface)  serves  as  a  transition 
region  from  plane  stress  to  plane  str.airi,  while  the  remaining  70%  is  under 
nominally  plane  strain  conditions. 

Figure  4.7  presents  the  residual  zero  load  mid-plane  opening  obtained 
from  optical  interference  for  increasing  crack  length.  Note  here  that  for  a 
crack  size  of  10.3  mm  the  point  of  maximum  opening  (bump)  is 
approximately  2  mm  away  from  the  mid-plane  crack  tip,  and  when  the 
crack  length  reaches  19.1  mm,  this  bump  lies  4.4  mm  away  from  the  local 
crack  tip.  As  a  result,  an  increase  of  8.8  mm  in  crack  size  corresponds  to  an 
increase  of  2.4  mm  in  the  distance  between  the  crack  tip  and  the  point  of 
maximum  opening.  Thus  this  ’bump’  where  the  crack  surface  separation  is 
a  maximum  under  zero  load,  moves  in  the  direction  of  the  advancing  crack 
tip. 


4.1.2  CMOD  and  BFS  Measurements: 


Ar  Hiscussed  earlier,  CMOll  and  BFS  ineasuretiienls  were  also  employed 
to  determine  the  averaj^e  (^lohal)  closure  loads.  For  the  ('MOD  ^a^e 
rn(>asiirement,  a  clip  was  mounted  ;il  the  starter  notch  mouth  .and  tin; 

signal  from  the  gage  was  recorckul  as  a  function  of  .ap()ru!d  load.  Figure  F8 
shows  the  CMOD  gage  reading  versus  applied  load  for  an  average  crack 
length  of  16.8  mm.  The  left-most  graph  prr'sents  the  load  versus 
displacement  curve  obtained  from  the  clip  gage  at  the  crack  moiilh  during 
one  cycle  of  loading  (The  right-hand  curve  is  discussed  later).  Note  that  the 
displacement  varies  in  a  non-linear  fashion  up  to  a  certain  load,  followed  by 
larger  linear  regioti,  and  a  slight  non-linear  region  at  the  upper  end  of  the 
load  displacement  record.  As  discussed  earlier,  the  closure  load  may  l)e 
determined  from  this  data  as  the  point  where  the  load-displacement  relation 
first  becomes  linear.  The  small  non-linear  region  at  the  upper  end  of  the 
load  cycle  may  be  due  to  large  plastic  deformation  at  the  crack  tip  present 
at  the  high  load  levels  of  the  cycle.  Moreover,  Figure  4.8  indicates  that 
PMiVlA  exhibits  fairly  large  hysteresis.  P’or  the  present  study,  the  loading 
curve  is  used  to  determine  the  closure  loads. 

The  closure  load  is  determined  here  by  an  olfset  axis  method,  although 
there  are  various  other  techniques  described  in  the  literature  to  obtain  this 
closure  value.  In  this  offset  axis  technique,  one  least  sfiuares  line  is 
calculated  through  the  upper  linear  portion  of  the  loading  curve.  The 
deviation  of  the  actual  displacement  from  the  least  squares  line  is  plotted  in 
an  expanded  scale  (approximately  5X)  in  the  right-half  of  Figure  4.8.  The 
closure  load  is  then  determined  as  the  point  where  the  actual  data  deviates 
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from  the  least  squares  line,  and  is  referred  to  here  Kq^.  Here  Kqc  's  defined 
as  the  crack  opening  load  obtained  from  the  ('MOD  measurements,  and 
from  figure  l.S,  for  an  average  crack  size  of  16.8  mm,  K()(;  /  '(lualed 

(),:o. 

A  similar  procedure  was  used  to  obtain  the  crack  closure  load  from  the 
back  face  strain  measurements  as  seen  in  f'igure  4.9.  Here  the  readings  from 
a  1000  ohm  strain  gage  mounted  on  the  back  face  of  the  specimen  were 
analyzed  to  determine  a  closure  load.  This  measurement  gave  a  /  Hmax 
value  of  0.34,  where  Kqe  is  defined  as  the  crack  opening  load  determined 
from  the  BPS  data. 

Figure  4.10  compares  the  closure  loads  determined  by  interferometry, 
C.VIOD,  and  BFS  measurements  for  increasing  crack  length.  Here  the  crack 
opening  loads  increase  slightly  for  increasing  crack  length  up  to  a  point 
before  returning  to  a  lower  value  for  further  crack  extension.  Note  from 
Figure  4.10  that  the  interferometric  free  surface  opening  load  is  larger  than 
the  compliance  (CMOD  and  BFS)  opening  load.  Since  interferometric 
measurements  showed  that  the  crack  tip  in  the  specimen  interior  was 
separated  under  zero  load,  the  global  closure  loads  obtained  from  the  CMOD 
and  BFS  techniques  should  fall  between  zero  and  the  interferometric  free 
surface  opening  load. 


4.1.3  Sensitivity  Analysis: 

The  objective  of  this  section  is  to  describe  the  effects  of  specimen 
thickne.ss,  cyclic  and  stress  ratio  (R)  on  the  fatigue  crack  closure 


behavior  of  through-thickness  flaws  in  the  PMMA  specimens.  Several 
experiments  conducted  in  this  test  program  are  stimrnarized  in  d'abU'  1.1. 

4. 1.3.1  Thickness  Effect: 

Five  fatigue  crack  closure  experiments  were  conducted  to  determine  the 
etfect  of  specimen  thickness  on  crack  opening  behavior.  The  cyclic  for 

these  tests  was  ke;)t  a  constant  at  hOO  Kl*a-m*'",  and  the  stress  ratio  (11) 
was  kept  at  0.1,  while  the  specimen  thickness  ranged  from  5  mm  to  25  mrn. 
,\s  indicated  in  'I'able  1.1,  the  specimen  thicknesses  for  tests  PX-16,  P.V-13, 
PX-12,  PX-11,  and  [*X-10  were  5,  10,  13,  19,  and  25  mm  respectively.  Uecall 
that  the  PM.M.A  specimens  were  cut  from  a  single  25  mm  thick  sheet.  The 
19  mm  thick  test  specimens  were  prepared  by  removing  a  6  mm  layer  of 
material  from  one  side  of  the  25  mm  thick  members.  However,  the  thinner 
test  specimens  were  prepared  by  removing  an  equal  amount  of  material  from 
both  sides  of  the  25  mm  thick  sheet.  As  mentioned  earlier,  all  specimens 
were  annealed  after  machining  to  the  prescribed  thickness. 

Figure  4.11  shows  the  crack  tunnelling  behavior  for  various  thicknesses. 
Here  the  difference  between  the  free  surface  and  mid-plane  crack  lengths  is 
plotted  against  the  average  (compliance)  crack  size.  Note  from  Figure  1.11 
that  the  amount  of  crack  tunnelling  remains  relatively  constant  with  respect 
to  increasing  crack  size,  but  varies  significantly  with  changing  si)ec,iinen 
thickness.  For  example,  the  ditference  between  the  free  surface  and  the 
mid-plane  crack  dimensions  for  a  25  mm  thick  specimen  under  steady  state 
cycling  is  about  4.5  mm,  while  for  a  5  mm  thick  member  the  magnitude  of 
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crack  tunnelling  is  negligible. 

I'igure  1.12  <‘onipar('s  the  /.('ro  load  optical  interference  fringe  patterns 
for  tlie  2’)  mm  am!  a  mm  tliick  im-mlx'rs.  .\'ote  here  that,  although  a  major 
portion  dl'  both  crack  planes  remain  opened  iiiid<-r  zero  load  for  both 
members,  the  crack  tii)  in  the  mid-plane  remains  partially  closed  in  the  .5 
mm  thick  specimen.  Furthermore,  the  zero  load  fringe  patterns  for  the  o 
mm  thick  rtuunber  are  <pjite  ditferent  from  the  25  mm  thick  specimen.  The 
fringes  in  the  25  mm  thick  member  form  an  elliptical  pattern  in  the  middle 
portion  of  the  crack  plane  indicating  a  large  displacement  gradient  between 
the  mid-plane  and  the  free  surface  opening  behavior,  while  for  the  5  mm 
thick  member,  the  fringe  patterns  are  relatively  straight  across  the  crack 
plane  indicating  a  more  uniform  through-thickness  opening  behavior. 
Furthermore,  the  25  mm  thick  specimen  clearly  shows  the  closed  portions  of 
the  crack  plane  near  the  free  edges  under  zero  load,  while  a  major  portion  of 
the  free  surfaces  for  the  5  mm  thick  member  remain  open  under  zero  load. 

Figure  t.1.1  presents  these  zero  load  (residual)  displacement  fields 
quantitatively  for  various  specimen  thickness.  Here  the  mid-plane 
displacement  is  plotted  in  terms  of  fringe  order  as  a  function  of  distance 
from  the  crack  tip.  Note  in  Figure  4.13  that  the  magnitude  of  the  mid-plane 
opening  decreases  with  a  reduction  in  the  specimen  thickness.  Furthermore, 
the  zero  load  opening  profile  for  the  5  mm  member  reflects  the  closed  portion 
of  the  crack  at  the  mid-plane  location. 

Figure  4.M  presents  the  interferometric  opening  load  for  different 
locations  through  the  specimen  thickness.  Note  from  figure  4.i4  that  when 
the  absolute  distance  (x)  from  the  free  surface  is  considered,  the  width  of  the 
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transition  region  remains  invariant  with  the  specimen  total  thickness.  For 
example,  in  figure  f.lt,  the  transition  region  is  lonfineil  to  a  stri[) 
approxi//i.af('ly  I’-.'f  mm  from  c.ach  i'ri  <•  surface  regardless  of  the  s[>eclmen 
thickness.  This  width  of  i  tie  Ir.an  -ilion  region  is  consistent  with  the  fact 
that  for  a  .a  mm  thick  specitnen,  })ortions  of  the  mid-plane  crack  lip 
remain(.'d  closed  under  zero  load  indicating  that  the  entire  specimen 
thickness  served  as  the  transition  region  for  this  thin  member. 

Figure  1.15  presents  the  interferometric  free  surface  opening  load  as  a 
function  of  the  average  (compliance)  crack  length  for  the  various 
experiments  discu.ssed  here.  Note  here  that  tlu'  opening  loads  do  not  differ 
significantly  for  the  25  trim  and  If)  mm  thick  members,  however,  these  free 
surface  opening  values  decrease  with  further  reduction  in  the  specimen 
thickness.  This  decrease  in  crack  opening  load  for  a  'thin'  member  is  not 
consistent  with  the  plasticity  induced  closure  arguments.  Since  a  major 
portion  of  a  thin  member  is  under  plane  stress  conditions,  the  crack  opening 
loads  for  these  thin  members  should  be  higher  than  the  thick  specimens  due 
to  the  presence  of  higher  compressive  residual  stresses  in  the  large  plastic 
zones  of  the  thin  members.  Since  the  experimental  results  showed  that  the 
crack  opening  load  decreases  with  a  sufTicient  decrease  in  specimen 
thickness,  perhaps  some  other  factors  such  as  the  crack  surface  roughness, 
which  Increases  with  a  reduction  in  specirtum  thickness,  plays  an  important 
role  in  the  crack  opening  behavior  of  the  thin  members.  However,  this 
phenomenon  where  the  crack  opening  load  decreases  with  a  reduction  in 
specimen  thickness  is  not  clearly  understood  at  the  present  time. 


4. 1.3. 2  Cyclic  and  Stress  Ratio  (R)  Effect: 


'fho  objective  of  lliis  sect 'on  is  to  describe  the  etlecls  of  cyclic 
the  str(‘ss  rntio  on  tlu'  crack  openitej,  Ix-ltavior,  As  sliown  in  'I'abh'  1.1,  tlu.' 
cyclic  for  the  I!)  tnrti  thick  members  I’.X-ll,  I’X-Ol',  I’X-O.'i.  and  I’.X-Ol 

were  (iOO,  hhO,  82'),  and  9,'b')  Klhi-m'-'',  and  for  the  2')  mm  thick  members 
PX-15,  PX-10,  PX-(),5,  PX-()7,  and  PX-l  l  were  110,  (>00,  715,  715,  and  770 
IXPa-m^^"  respectivady,  .Altiiough  the  cyclic  ("or  tests  PX-05  and 

PX-07  were  kept  a  constant  at  715  KPa,-m''",  the  If  values  for  thes('  tests 
were  0.1  and  O.,!  respect! velja 

Figure  4.16  summarizes  the  interferometric  fr('e  surface  opening  loads 
for  the  variou,s  experiments  with  19  mm  thick  members,  and  Figure  1.17 
presents  the  results  for  the  25  mm  thick  specimens.  Note  from  these  figures 
that  the  crack  opening  loads  for  the  intermediate  cyclic  experiments 

are  somewhat  higher  than  the  experiments  conducted  under  relatively  Mow’ 
or  ’high’  levels.  For  example  in  Figure  1.17,  the  crack  opening  loads 

observed  under  steady  state  cycling  of  600  KPa-m‘'^^  are  higher  than 

the  experiments  conducted  under  cyclic  of  410  and  770  KPa-rn''^’. 

Although,  it  is  expected  that  the  crack  opening  loads  should  decrease  with  a 
redtiction  in  the  cyclic  level  from  plasticity  argtiments,  another 

explanation  is  in  order  for  the  reduction  in  crack  opening  load  at  extremely 
high  cyclic  I’erhaps  the  ditference  in  the  crack  closure  mechanisms 

with  varying  cyclic  K  ranges  combined  with  the  rough  crack  surfaces 
observed  under  high  K  levels,  contribute  to  thi.s  crack  opening  behavior.  For 
example,  plasticity  may  dominate  the  crack  opening  behavior  at  ’low’  and 
’intermediate’  cyclic  ranges,  while  crack  surface  roughness  combined 


with  crack  craze  zone  bundling  have  significant  effect  on  the  closure  values 
under  'high’  cyclic  stress  levels.  'I'he  fact  that  the  crack  opening  loads 
decre.ases  with  an  iticrease  in  cyidic  Uwel.  is  consistent  \slth  iIh'  resui’s 

in  the  [)revious  s(_'clion,  wlu're  the  crack  o|)<'ning  loads  decrease  with  a 
suHicient  decrease  in  specimen  thickness  du(>  to  the  crack  surl'ace  rougliness 
arguments.  Sincc^  the  free  surface  crack  separation  lo.ad  attained  a 
ma.ximum  level  for  a  specimen  subjected  to  an  intermediate  cyclic  the 

cr.ack  tunnelling  magnitude  under  this  level  is  significantly  higher  than 

the  specimens  subjected  to  either  ’high’  or  ’low’  cyclic  stress  values  as  seen 
from  Figure  4.18. 

Figure  4.19  presents  the  interference  opening  load  for  various  distances 
from  the  specimen’s  free  surface  for  the  25  nrun  thick  members.  .\ote  here 
that,  although  the  free  surface  opening  loads  vary  somewhat  with  the  cyclic 
^max  levels,  the  width  of  the  transition  region  (near  the  free  surface)  does 
not  vary  significantly  with  the  K^ax-  Fo*"  example,  the  closure  phenomenon 
is  confined  to  a  region  about  2-3  mm  from  the  free  edges. 

Figures  4.20  and  4.21  compare  the  free  surface  opening  loads  with  the 
opening  values  obtained  from  CMOD  and  BFS  readings  for  tv. -^s  PX-0.5  and 
PX-07  respectively.  Recall  that  although,  the  cyclic  levels  for  these 

two  experiments  were  set  at  715  KPa-m*'*^^,  the  stress  ratio  (R)  for  these 
tests  were  kept  a  constant  at  O.l  and  0.3  for  the  tests  PX-05  and  PX-07 
respectively,  'rhus  for  Test  PX-05  was  set  at  72  KPa-m'^'  while  it 

equaled  215  KPa-m*^^  for  Test  PX-07. 

Note  from  figures  4.20  and  4.21  that,  although  the  crack  o[)oning  levels 
for  these  two  experiments  do  not  vary  significantly,  the  crack  o[)ening  loads 
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for  Test  PX-07  (R  —  0.3)  are  less  than  or  near  the  ryrlic  so  that  th(.‘ 

crack  is  almost  always  open  dnrin*'  the  l)aselin('  cycling.  'I'his  f.ict  is 
rellected  in  the  <)l)S(‘rvat  ion  that  the-  R  0.1  cr.’uk  ”ro\‘.  t.h  r;ilc  (l..')7  \  10  ' 
mrn/cych?)  w;is  less  lhan  the  It  0.3  cr.ai'k  'growth  r;it('  (l.i)t  \  10  * 
inin/cy<‘le),  althougli  the  liaseline  cyclic  AK  for  the  R  0.1  lest  was  higher 
(643  KPa-ni'''*)  than  the  R  =  0.3  e.xperiinent  (AK  —  .300  KPa-m''").  Thus 
the  effective  cyclic  AK  for  these  two  experiments  were  nearly  etpial  in 
magnitude,  (free  surface  AK  -  lid  KPa-m'^“)  althoirgh  the  nomin.al  applied 
cyclic  AK  for  tests  PX-5  and  PX-7  were  kept  at  643  and  .300  KPa-m'  " 
respectiv'ely. 


4.2  THROUGH-THICKNESS  FLAWS  IN  PC: 

The  objective  of  this  section  is  to  describe  the  three-dimensional  aspects 
of  fatigue  crack  closure  observed  in  the  polycarbonate  members.  As  before, 
the  point-wise  measure  of  crack  opening  load  obtained  via  optical 
interferometry  is  compared  with  the  bulk  measurements  from  CMOD  and 
BPS  reading.  A  total  of  six  through-thickness  flaw  experiments  were 
conducted  with  the  polycarbonate  material  as  summarized  in  Table  4.2. 
Figure  4.22  compares  the  crack  tunnelling  behavior  among  the  25,  19,  and  10 
mm  thick  specimens.  Note  here  that  as  with  the  through-thickness  Haws  in 
the  PMMv\  members,  the  magnitude  of  crack  tunnelling  decreases  rapidly 
with  a  reduction  in  specimen  thickness. 

Figure  4.23  shows  the  crack  opening  perimeter  as  a  function  of  applied 
load.  Again  the  area  enclosed  between  the  crack  front  and  the  boundary 
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line  for  K  /  =  0  represents  the  initial  (residual)  opened  section  of  the 

crack  plane.  .\ote  in  Figure  l.'id  that  unlike  the  I’MMA  experiini.uits,  both 
the  free  surface  and  niid-{)lane  locations  of  this  'JT)  nun  thick  P('  ineiuber 
remain  opened  under  zero  load,  howev<'r,  an  area  of  the  crack  [)erim('t (t  in 
the  interior  of  the  specimen  remained  closed  under  zero  load.  'I'he  residual 
free  surface  opening  may  suggest  that  the  closure  phenomenon  is  confined  to 
a  region  near  each  free  surface,  but  not  at  the  free  surface  in  these  thick  PC 
test  members.  This  crack  opening  behavior,  wlusre  the  crack  opens  last  at  a 
point  interior  to  the  free  surfaces,  was  also  observed  in  surfacm  flaw 
experiments  subjected  to  cyclic  bending  in  reference  |48j.  Figure  4.24  shows 
the  crack  opening  perimeter  as  a  function  of  applied  load  for  a  19  mm  thick 
member.  The  crack  opening  perimeter  for  this  19  mm  thick  specimen  follows 
the  same  pattern  observed  in  the  PM.VIA  experiments,  where  the  maximum 
closure  occurred  at  the  free  surface  locations. 

Recall  from  above  that  the  test  member  PC-3  was  prepared  by 
removing  a  3  mm  layer  of  material  from  both  sides  of  a  25  mm  thick  sheet, 
while  test  pieces  PC-2  and  PC-5  were  constructed  by  removing  a  6  mm  layer 
from  one  side  of  the  25  mm  thick  member.  Although  these  specimens  were 
cut  differently,  the  crack  opening  profiles  observed  in  these  three  tests  were 


similar. 


Figures  4.25  and  4.26  present  the  digitized  crack  opening  profiles 
obtained  at  the  rnid-plane  and  free  surface  locations  respectively  for  a  19 
mm  thick  member.  Here  the  magnitude  of  crack  opening  in  the  specimen's 
interior  is  significantly  larger  than  the  free  surface  opening,  consistent  with 
the  PMMA  experiments  described  earlier.  This  difference  in  crack  opening 
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between  the  free  surface  and  the  mid-plane  is  presumably  due  to  the 
variation  of  plasticity  present  along  the  crack  perimeter.  I'igures  1.27  and 
1.28  show  the  crack  opening  l)(diavior  for  a  25  mm  thick  member  recorded  at 
the  mid-plane  and  near-free  surface  locations  respc-ctivcdy.  Rc-call  that  in 
the  25  mm  thick  members,  the  free  surface  remained  open  under  zero  load. 
Figure  ‘1.28  shows  the  crack  opening  behavior  at  a  depth  x/B  ^  0.08,  w}}ere 
X  is  the  distance  from  the  free  surface,  and  B  is  the  specimen  thickness. 

Figures  4.29  and  4.30  present  the  interferometric  crack  tip  opening  load 
as  a  function  of  distance  from  the  free  surface  for  the  25  mm  and  19  mm 
thick  specimens.  Here  the  distance  from  the  free  surface  (x)  is  normalized 
with  the  specimen  thickness  (B)  and  the  opening  load  is  normalized  with  the 
cyclic  maximum  load.  Note  in  these  figures  that  as  before,  the  crack  tip 
separation  load  decreases  rapidly  as  one  moves  into  the  specimen  thickness. 

Figure  4.31  compares  the  interferometric  free  surface  opening  loads  with 
the  bulk  measurements  obtained  from  the  CMOD  and  BFS  reading.  Recall 
that  the  bulk  opening  loads  from  the  CMOD  and  BFS  reading  are 
determined  as  the  point  where  the  load-displacement  curve  first  become 
linear  during  the  initial  opening  phase  of  a  load  cycle.  Note  from  figure  4.31 
that  these  bulk  measiirements  are  again  generally  lower  than  the 
interferometric  free  surface  opening  load,  although  the  difference  between 
these  loads  are  considerably  less  than  the  results  obtained  from  the  PMMA 
experiments,  where  the  bulk  opening  loads  were  found  to  be  significantly  less 
than  the  free  surface  opening  load. 


P’igure  4.32  compares  the  zero  load  mid-plane  opening  profile  for  three 
different  specimen  thicknesses.  Note  here  that  the  magnitude  of  the  residual 
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crack  surface  separation  again  decreases  with  a  reduction  in  specimen 
thickness,  conslstcMit  with  the  behavior  exhibited  with  the  f^MMA  tiirough- 
ilaw  ('X[)('rinieiit.s. 

figure  presents  the  interferometric  free  surface  opening  loarls  for 

tile  various  polycarbonate  members  reported  here.  Here  the  opening  loads 
from  the  25  and  19  mm  specimens  form  a  scatter  band  between  50  and  115 
KPa-m*^^.  The  baseline  cycling  w’as  at  of  330  KPa-m''^“  with  R  equal 

to  0.1  for  all  cases. 


4.3  CORNER  CRACKED  HOLES: 

Table  4.3  presents  the  test  matrix  for  the  corner  cracked  hole 
experimental  program.  The  cyclic  maximum  remote  tensile  stress  for  tests 
CTB-1  and  CTB-4  was  kept  at  3150  KPa,  and  for  tests  CTB-3  and  CTB-5 
was  at  4200  KPa.  Following  the  tensile  portion  of  the  experiment,  the 
maximum  cyclic  bending  stress  for  tests  CTB-1  and  CTB-2  was  maintained 
at  7360  KPa,  and  for  tests  CTB-3  and  CTB-5  was  kept  at  5520  KPa.  The 
stress  ratio  (R)  for  all  tests  was  maintained  at  0.1. 


4.3.1  Cyclic  Tension  Results: 

Optical  interference  was  employed  to  obtain  three-dimensional  crack 
surface  displacement  profiles  as  a  function  of  applied  load.  Figure  4.34 
shows  the  fringe  patterns  obtained  during  steady  state  cycling  for  an 
increasing  remote  tensile  stress  when  the  crack  dimensions  a  7.7  mm  and 
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c  =  4.6  nim.  Figure  4.54(a)  shows  the  interference  pattern  at  zero  load,  and 
indicates  a  residual  opening  displac(unerit  field.  Recall  that  ex{)eriiiients 
conducted  with  through-flaws  in  hotfi  f'.\I\f.\  and  FC"  edge  cr.icknMl  mi-nihers 
showed  a  similar  residual  patterti  uiuler  zero  load.  In  I'igui-e  1.51(:i)  tlu' 
fringes  in  the  middle  portion  of  the  <-rack  plane  indicate  that  particular 
section  of  the  crack  plane  is  open  under  zero  load,  while  the  crack  surfaces 
at  0  =  0°  and  90°  (the  locations  where  the  flaw  intersects  the  free  surface 
and  the  hoh;  bore)  rcmiairi  closed.  As  the  remote  tensile  load  is  increasi'd, 
the  0-order  light  fringe  reaches  the  0°  and  90°  locations  at  a  stress  level  of 
331  KPa.  Although  the  a  and  c  crack  locations  open  at  the  same  remote 
stress  level,  the  stress  intensity  value  (K)  at  these  two  points  are  different. 
Here  the  stress  intensity  factor  for  a  particular  location  is  computed  from 
the  Newman  and  Raju  solution  [57j  employing  the  nominally  applied  remote 
tensile  stress. 

Figure  4.35  shows  the  opened  crack  perimeter  as  a  function  of  applied 
remote  tensile  stress.  Recall  that  Figure  4.35  is  obtained  by  overlapping  the 
opened  areas  of  the  crack  for  increasing  load  for  a  particular  crack 
geometry.  Thus  the  area  enclosed  between  the  original  crack  front  and  the 
boundary  line  for  for  a  particular  stress  rr  represents  the  opened  section  of 
the  crack  plane  under  the  respective  applied  remote  stress.  Note  here  that 
the  crack  front  first  starts  to  open  along  the  middle  portion  of  the  corner 
flaw,  and  then  progresses  outward  to  the  hole  and  the  free  surface  locations. 

The  fact  that  a  larger  load  was  required  to  open  the  crack  faces  at  the 
free  surfaces  than  at  interior  locations,  was  also  observed  for  the  through- 
crack  experiments  reported  in  earlier  sections.  Those  experiments  indicated 
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that  the  crack  opens  first  at  the  specimen’s  mid-plane  followed  by  opening 
at  the  free  surfaces  at  a  signiiicanlly  higher  load.  This  difference  in  crack 


ofauiing  is  ()reMiiiial)ly  due  l(i  the  variation  in  plastic  /one  size  along  tin- 
cr.ick  front,  as  c<)iiditions  change  from  plane  str(‘ss  to  plane  strain. 

Crack  opening  displacements  were  examined  (piantitatively  by  tlu; 
digitized  measurements  of  the  interference  fringe  patterns  of  Figure  4. .31. 
Here  the  displacements  are  presented  in  both  fringe  order  (n)  and 
corresponding  metric  (/on)  units.  Figure  4.36  presents  the  crack  surface 
separation  under  a  given  remote  stress  (cr)  for  varying  parametric  angle 
In  this  figure,  the  local  original  crack  tip  serves  as  the  plot  origin.  Here  the 
displacement  profiles  are  obtained  perpendicular  to  the  local  crack  front  at 
various  points  along  the  crack  perimeter  defined  by  the  elliptic  angle  o.  For 
example,  the  displacement  curve  for  </>  =  30°  represents  the  crack  face 
separation  along  a  line  perpendicular  to  the  crack  front  at  a  point  defined 
by  the  elliptic  angle  d’  =  30°.  Note  from  Figure  4.36  that  the  magnitude  of 
the  displacement  field  tends  to  increase  as  one  moves  away  from  either  the 
specimen’s  free  surface  or  the  hole  location.  The  displacement  curves  for  d' 
=  0°  and  90°  in  Figure  4.36  indicate  that  those  portions  of  the  corner  flaw 
are  closed  under  the  applied  remote  tensile  stress  of  263  KPa. 

Figure  4.37  presents  the  crack  opening  profile  measured  at  the  free 
surface  -  -  0°).  Fach  curve  here  represents  a  different  applied  remote 
stress,  and  gives  the  total  crack  separation  (2V)  between  the  two  crack 
surfaces  as  a  function  of  distance  from  the  original  crack  tip.  Recall  that, 
the  stress  which  causes  the  displacement  curve  to  pass  through  the  origin, 
giving  complete  crack  separation,  is  referred  here  as  the  interferometric 


opening  load  for  the  particular  crack  location.  In  Figure  4.37  the  crack  tips 
separate  completely  at  the  specitnen’s  free  surface  under  an  applied  local  K 
of  t)7  KP,i-in'  I’igure  1.3S  presetils  the  corresponding  crack  openiing 
protile  iii<'aTure<l  along  the  h.ole  i)or<‘.  'I'lie  opening  K  value  determined  from 
Figure  l.liN  ecpials  71  Kl'a-m^'“.  Thus  portions  of  the  interior  crack  plane 
are  openetl  under  zero  load,  wliile  the  free  surface  and  hole  locations  require 
opening  loads  of  67  and  7l  KPa-in*'^"  respectively.  It  is  important  to  note 
that  file  !(;cal  K  values  are  computed  froni  the  .\evvman  and  Raju  solution 
r)7j,  and  these  de[)end  on  '•  as  well  as  load.  Thus  for  a  particular  bending 
.stress,  K  vari<‘s  along  the  crack  perimeter  with  the  elliptic  angle 

Figure  4.39  presents  the  interferometric  crack  tip  separation  load 
determined  at  different  parametric  angles  for  test  CTB-1.  Here  the  opening 
loads  are  normalized  with  the  cyclic  maximum  remote  load.  Note  that  the 
opening  load  decreases  rapidly  as  one  moves  away  from  either  the  free 
surface  or  the  hole  location.  This  decrease  in  opening  load  with  increasing 
distance  from  the  free  surfaces  was  also  observed  with  the  through-thickness 
flaw  experiments  reported  earlier. 

4.3.2  Cyclic  Bend  Results: 

Figure  4.10  presents  a  typical  set  of  fringe  pattern  photographed  during 
the  opening  phase  of  a  bending  load  cycle.  Figure  4.41  shows  the 
corresponding  crack  opening  perimeter  as  a  function  of  applied  remote 
bending  stress.  Note  in  Figure  4.41  that,  unlike  the  opening  behavior  under 
remote  tension,  the  free  surface  opens  earlier  than  the  hole  bore  location. 


For  example,  the  crack  surfaces  at  <!>  ■=--  0°  separate  under  a  remote  bending 
stress  of  8(52  Kf’a,  while  tlu;  crack  surface's  along  the  hole  bore  require  an 
applied  stress  of  1731  KFa  to  separate.  'I'liis  dilference  in  opening  loads  is 
caused  by  the  sharp  stress  gradient  pres(uit  along  the  hole  direction. 

Figure  4.42  presents  the  crack  opening  profiles  for  various  parametric 
angles  o  under  a  given  remote  bending  stress  of  862  KPa.  Note  here  that 
the  displacement  magnitude  again  increases  as  one  moves  away  from  either 
the  free  surface  or  hole  locations.  Recall  that  similar  opening  behavior  was 
observed  earlier  for  the  remote  tension  phase  of  the  experiment.  Here  the 
free  surface  displacement  field  at  ■;>  —  is  significantly  larger  than  the  90° 
profile  due  to  the  fact  that,  under  a  remote  bending  stress  of  862  KPa,  the 
K  level  experienced  at  the  free  surface  location  is  significantly  higher  than 
the  K  value  at  the  90°  crack  tip  location.  For  the  crack  dimensions  under 
consideration,  a  remote  bending  stress  of  862  KPa  corresponds  to  K  levels  of 
123  and  42  KPa-m’^^  at  the  free  surface  and  the  hole  bore  locations 
respectively. 

Figure  4.43  presents  the  crack  opening  profile  for  increasing  remote 
bending  stress  at  the  free  surface  location  {<!>  —  0°).  Note  here  that  an 
applied  bending  stress  of  862  KPa  (K  =  123  KPa-m*'^^)  is  needed  for  the 
crack  opening  curve  to  pass  through  the  origin  (crack  tip)  in  order  to  give 
complete  crack  tip  separation.  Similar  crack  opening  profile  plots  obtained 
along  the  hole  bore  (c'>  —  90°)  are  shown  in  Figure  4.44.  A  remote  bending 
stress  of  1734  KPa  (K  =  84  KPa-m'^^)  is  required  here  to  open  the  crack  tip 
completely  along  the  hole  bore. 


Figure  4.45  presents  the  remote  bending  opening  stress  (/r^)  as  a 
functit)n  of  parametric  angle  'I’he  ojxming  stress  is  normalized  here  with 


ri's()ect  to  the  m.aximum  cyclic  bendiiig  stress  .\oLe  in  these* 

figures  that  as  tin*  crack  grows,  the*  '’‘o/'ani<  value  is  reduced,  and  a  larger 
portion  of  the  crack  front  remains  propped  open  under  zero  load.  For 
example,  the  residual  displacement  field  for  a  crack  with  a  =  7.7  and  c  = 
4.6  mm  is  smaller  than  the  residual  field  for  a  crack  with  a  =  7.7  and  c  = 
8.8  mm. 


4.4  SURFACE  FLAWS  UNDER  CYCLIC  BENDING: 

The  objective  of  this  section  is  to  describe  experiments  directed  at 
determining  three-dimensional  aspects  of  fatigue  crack  closure  for  a  surface 
flaw  in  a  rectangular  plate  subjected  to  pure  bending.  As  before,  optical 
interference  technique  is  employed  here  with  the  transparent  PMMA 
specimens  to  determine  the  point-wise  measurement  of  fatigue  crack  closure. 
These  point-wise  measurements  are  then  compared  with  an  average  (bulk) 
crack  opening  load  obtained  via  CMOD  gage  measurements.  The  test 
matrix  is  shown  in  Table  4.4. 


4.4.1  Interferometric  Measurements: 

Figure  4.46  presents  a  set  of  interference  fringe  patterns  photographed 
for  increasing  load  for  tests  SB-2  with  crack  dimensions  a  =  5.8  mm  and  2c 
=  13.5  mm.  Here  cr  is  the  maximum  remote  3tre.s3  computed  by  the 


a-' 


SI 


V" 


ft 


V 

if 


•-  .V 


"•S' 

i 


nV 


v^y 


r-.v 


3'W 


■'•‘a 


. 


61 


standard  flexural  stress  formula.  The  fringe  pattern  at  zero  load  again 
indicates  a  residual  crack  opening  displacement  field. 

Note  from  Figure  4.46  that  under  zero  load,  fringes  are?  pri'seiU  in  the 
middle  portion  of  the  surface  Haw,  indicatitig  that  section  of  (  he  crack  plain- 
is  open,  while  the  free  surface  crack  tip  locations  require  a  stres.s  level  of 
2482  KPa  to  open.  Thus  the  effective  cyclic  stress  varies  along  the  crack 
front.  Figure  4.47  shows  the  opened  crack  perimeter  as  a  function  of  applied 
remote  stress.  /Vs  before.  Figure  4.47  is  obtained  by  overlapping  the  opened 
areas  of  the  crack  from  Figure  4.46  for  increasing  applied  load. 

Figures  4.48  and  4.49  present  the  crack  opening  profiles  at  various 
points  along  the  crack  front.  Here  the  crack  surface  displacement  is 
presented  in  fringe  order  units  as  a  function  of  distance  from  the  original 
crack  tip.  Figure  4.48  shows  the  displacement  profiles  along  slices  made  at 
different  parametric  angles  (0)  under  an  applied  load  of  20%  of  the 
maximum  remote  cyclic  stress.  These  slices  were  oriented  in  the  radial 

direction.  For  example,  the  curve  for  <■>  =  17°  in  Figures  4.48  and  4.49 

represents  the  crack  surface  displacement  profile  along  a  radial  line 

connecting  the  crack  origin  (X  =  Y  =  0.0  in  Figure  2.5)  and  the  point  on 
the  crack  front  where  equals  17°.  Recall  that  for  the  corner  cracked  hole 
specimens,  the  crack  surface  displacements  were  measured  perpendicular  to 
the  local  crack  front.  Note  in  these  figures  that  the  crack  surface 

displacements  tend  to  increase  as  one  moves  away  from  the  specimen’s  free 
edge  (</>  =  0°).  The  displacement  curve  for  (j)  =  0°  indicates  that  the  region 
near  the  crack  tip  is  closed  under  an  applied  load  of  20%  of  the  maximum 
cyclic  remote  stress.  Figure  4.49  shows  the  opening  profiles  for  a  load  level 
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which  is  28%  of  the  maximum  remote  cyclic  load,  and  indicates  that  the 
entire  crack  is  opened  under  this  stress  level. 

Figure  ■}.'>()  shows  t  he  iiilerferomc'lric  crack  sc'p.arat  ion  lo;id  (leterinined 
at  differemt  parametric  an^hjs  for  lest  SB-2.  Here  the  local  opening  stress 
(rr^)  is  normalized  with  respect  to  the  maximum  cyclic  stress 
observed  earlier  for  through-thickness  and  corner  flaw  experiments,  Figure 
4.50  shows  that  the  opening  stress  decreases  rapidly  as  one  moves  away  from 
the  free  surface  location. 

Figure  4.51  presents  the  free  surface  interferometric  opening  stresses  (at 
<:>  =  0^)  obtained  at  various  crack  geometries  for  all  tests  reported  here. 
Recall  that  the  interior  ’a’  location  was  always  opened  under  zero  load, 
while  the  free  surface  required  a  significant  portion  of  the  load  cycle  to  open 
at  the  tips.  Thus  an  effective  cyclic  stress  level  can  be  computed  at  each 
point  where  the  ’a’  and  ’c’  dimensions  are  defined.  The  effective  stress  level 
at  a  given  point  is  defined  here  as  the  difference  between  the  maximum 
applied  cyclic  stress  and  the  stress  required  for  crack  tip  separation  at  the 
point  under  consideration. 

Note  in  Figure  4.51  that  the  free  surface  opening  loads  for  the 
intermediate  cyclic  experiments  are  somewhat  higher  than  the  opening 

loads  observed  under  relatively  ’low’  or  ’high’  cyclic  For  example,  in 

F'igure  4.51  the  crack  opening  loads  observed  under  cyclic  of  8790  KPa 
are  higher  than  the  opening  loads  obtained  from  experiments  conducted  at 
10,340  or  7240  KPa.  Similar  observations  were  made  with  the  through- 
thickness  flaw  experiments  where  the  crack  opening  loads  for  an 
intermediate  cyclic  were  higher  than  the  opening  loads  recorded  for 


either  ’low’  or  ’high’  cyclic  K^ax  experiments. 

[■'igure  l.o'i  shows  the  fatigue  crack  growth  rales  da/d.N  and  dc/<lf'' 
plotted  versus  the  corrc'sponding  nominal  applied  AK  for  the  five  cyclic  h('nd 
experiments  rf'ported  here.  Here  the  cyn-lic  st  ress  intensity  factors  n  re 
computed  from  the  Newman  and  Raju  solution  [59j  employing  the  iiomitially 
applied  cyclic  load.  A  least  square  fit  through  the  through-thickness  crack 
growth  rate  data  (da/dN  versus  nominal  applied  AK)  is  also  presented  in 
Figure  1.52  for  comparison  purposes.  Note  here  that  the  through  crack 
growth  behavior  lies  close  to  the  surface  Haw  growth  rates.  Figure  1.53  plots 
the  surface  crack  growth  rates  as  a  function  of  the  effective  AK  at  the 
respective  points  ’a’  and  ’c’.  Here  the  ’effective’  AK  is  computed  with  the 
’closure  corrected’  cyclic  load  which  results  when  the  crack  faces  are  always 
separated  at  the  points  of  interest. 

Note  from  Figure  4.52  that  da/dN  and  dc/dN  show  a  systematic 
variation  when  plotted  against  the  nominal  applied  AK.  For  example  in 
Figure  4.52,  a  least  square  power  law  (Paris)  fit  through  the  da/dN  versus 
applied  AK  data  shows  significant  deviation  from  a  similar  curve  fit  with 
the  dc/dN  quantities.  However,  when  the  closure  corrected  effective  AK 
value  is  used  in  place  of  the  nominal  applied  AK,  as  in  Figure  4.53,  the 
crack  growth  rates  da/dN  and  dc/dN  indicate  little  difference. 

4.4.2  CMOD  Measurements: 

As  discussed  earlier,  CMOD  readings  were  also  obtained  in  addition  to 
the  interferometric  measurements  during  the  crack  propagation  period.  The 


bulk  (average)  crack  opening  loads  were  determined  from  the  load- 
displacement  curves  obtained  from  a  clip  gage  mounted  at  the  crack  mouth. 

J’igure  1.5  1  conipares  the  crack  optuiing  loads  obtained  from  (d.MOD 
readings  with  the  interferometric  fret;  surface  opening  values.  Here  the 
opening  loads  are  plotted  versus  the  crack  aspect  ratios  (a/c)  to  observe  the 
effect  of  surface  flaw  shape  on  the  comparison  between  the  local  (free 
surface)  and  global  (CMOD)  opening  values.  The  dotted  lines  represent  the 
minimum  and  maximum  cyclic  stress  limits.  Hecall  that  at  the  specimen'.s 
depth  location  the  crack  tip  was  open  under  zero  load.  Note  from*  figure 
1.51  that  the  interferometric  free  surface  opening  load  generally  lies  above 
the  bulk  (average)  measurements  obtained  at  the  crack  mouth,  although  the 
difference  between  these  two  types  of  opening  loads  is  not  significant. 
Similar  comparisons  for  through-thickness  flaws  in  PMMA  members  using 
CMOD,  BFS,  and  optical  interference  measurements  showed  that  the  global 
crack  opening  loads  (from  CMOD  and  BFS)  were  significantly  lower  than  the 
interferometric  free  surface  opening  load.  Thus  the  bulk  crack  opening  load 
is  lower  than  the  free  surface  opening  value  for  through-thickness  flaws, 
while  for  surface  flaws,  these  two  quantities  do  not  vary  significantly. 

For  a  through-thickness  flaw,  the  entire  crack  front  is  equidistant  from 
the  CMOD  measurement  location;  thus  the  global  opening  load  obtained 
from  the  CMOD  reflects  an  average  of  the  mid-plane  and  free  surface 
opening  values.  On  the  other  hand  for  the  surface  (law,  the  C\fOD 
measurement  location  is  closer  to  the  free  surface  crack  tip  than  for  the 
maximum  crack  depth  point.  Thus  the  opening  load  determined  from 
CMOD  readings  should  be  significantly  influenced  by  the  free  surface  crack 
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tip  behavior,  and  may  not  reflect  the  true  global  (average)  quantity  for  th( 
entire  surface  flaw. 
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Tablf  4.1;  Tost  matrix  for  the  through-thickness  flaw  experiments  with 
P.MMA. 
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19 
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a  KPa 
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R 

Thickness  (mm) 

SB-1 
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25 

SB-2 
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0.1 

25 

SB-3 

10,340 

0.1 

25 

SB-6 

10,340 

0.1 

25 

SB-7 

7240 

0.1 

25 
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igure  4.4:  Mid-plane  crack  opening  profile  for  increasing  load.  The  crack 
opening  displacement  is  expressed  here  in  fringe  order  units. 
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Figure  1.10:  Comparison  of  crack  opening  loads  obtained  from  inter¬ 
ferometry,  CMOD,  and  BFS  measurements  for  various  crack 
lengths. 
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Figure  l.ll:  Comparison  of  the  magnitude  of  crack  tunnelling  with  chang¬ 
ing  specimen  thickness. 
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^'igiire  1.18:  C  ompari.'jori  of  crack  tunnelling  with  varying  cyclic  K 
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l'’lgiire  '1.21:  Crack  opening  load.^  obtained  from  various  mea.surcments  for 
increasing  crack  length  for  Test  PX-7,  R  =  0.3. 
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Figure  4.24:  Crack  opening  perimeters  as  a  function  of  applied  load  for  the 
19  rnrn  thick  PC  specimen. 
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Figure  1.27:  Mid-plane  crack  opening  profile  for  a  25  mm  thick  PC  speci¬ 
men  under  steady  state  cycling. 
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Figure  4.31:  Comparison  of  opening  loads  obtained  from  optical  inter¬ 
ferometry,  CMOD,  and  back  face  strain  measurements  for  Test 
PC-1  tested  with  R  =  0.1. 
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Fi<^urr  1.33:  Comparison  of  interferometric  free  surface  opening  load  fo 
varying  specimen  thickness. 
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Figure  4.:57;  Crack  opening  profiles  for  varying  applied  load  at  the  free  sur 
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Figure  -1.39;  Interferometric  crack  lip  opening  load  as  a  function  of  angle  <t> 
for  Test  CTB-1. 
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[’is^iir*'  1.12:  Crack  opening  profile  under  a  remote  bending  stress  of  86 
Kl’a  for  varying  parametric  angle  0. 
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Figure  4. HI:  Interferometric  free  surface  crack  opening  stresses  as  a  func 
tion  of  dimensionless  crack  size  for  five  cyclic  bend  expert 
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P"ig>ire  1.54:  Comparison  of  interferometric  free  surface  opening  load  with 
the  bulk  measurement  via  crack  mouth  opening  displacement 
for  tests  SB-1  and  SB-2.  The  maximum  remote  cyclic  bending 
stress  was  kept  at  8790  KPa. 
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ANALYTICAL  RESULTS 

Tlic  I'oinpiit (‘r  ;il^<jr'il  lim  (ic^scrihcd  in  K'linplcr  15  \v;is  n-cil  to  mcili  l  liic 
11  (iitl'Tcnt  surfaco  flaw  experiments  which  are  lislo<l  in  Table  a.l.  The 
calculation  was  performed  on  a  Cyber  205  and  running  limes  variiui  from 
aj)proximately  1000  to  1500  seconds.  The  running  lime  is  de[)endeni  upon 
the  applied  loads  and  the  amount  of  total  crack  growth.  'I'he  experiments 
include  five  bending  tests  by  Ray  [oOj  at  three  different  load  levads,  four 
bending  tests  by  Troha  ;62j,  and  two  tension  tests  by  Pope  ;58i.  .\ll  tests 
were  performed  on  polymethylmethacrylate  (PMXLV)  specimens.  Sections  5.1 
through  5.5  of  this  chapter  will  compare  some  of  the  experiments  with  their 
predictcfl  results.  Section  oA  will  examine  the  sensitivity  of  the  closure 
model  to  several  diiferent  parameters.  Section  4.1  of  reference  j58j  reviews 
the  mal('ri.al  [)roperlies  used  by  the  model,  which  will  not  be  gi\(m  in  this 
report.  Only  four  modfds  of  tlie  surface  Haw  experiments  are  summarized  in 
this  report,  .\ddition.al  rc'sults  for  the  other  experiments  arc  givcm  in 
reference  5.S  . 


6.1  TESTING  BY  RAY: 


As  slioun  ill  table'  .j.l,  l?ay  laO]  coruiuoted  five  bending  tests  at  three 
(lillerent  load  le\<  ls.  Optical  iiiterrero’net  ry  was  used  to  dete'riuiiie  when 
l■losurl■  occurred.  .\il  evperinients  wer*'  l•<>nducted  at  an  R  ratio  of  0.1. 

riu'  current  numerical  model  pre'dicts  significant  closure  at  th<'  interior 
of  tlie  surface  cracks,  which  wa.s  not  observed  in  any  of  the  exfieriments  with 
PM.\!.\.  Figure  5.1  shows  the  [iredicted  interior  closure  for  test  SB-1.  The 
iiioilel  showed  good  corre'latiou  with  the  experimt-ntal  results  for  closure  at 
the  fri'f'  surfaci'  for  tests  SB-l,  Slt-'i,  and  .SB-7  (figure  5.2  shows  the 
pr(  (iict«'d  free  surface  closure  for  test  SB-l).  'I'liese  tests  were  performed  at 
maximum  outi'r  iilier  bending  stresses  of  8790  and  7210  Kpa.  'Pests  SIT3 
and  SB-6  were  conducted  at  loaris  of  10,310  Kpa  and  exhibited  less  closure 
than  the  other  tests.  'Phe  free  surface  closure  predictions  by  the  model  were 
considerably  higher  than  the  experimentally  observed  values  for  tests  SB-3 
and  SB-6  (figure  5.3  shows  the  experimental  and  predicted  closure  for  test 
SB-3)  Because  of  the  interior  closure  calculated  by  the  model,  the  growth  in 
the  a  dimension  was  reduced,  and  the  aspect  ratio  (a/c)  of  the  modeled 
crack  became  less  than  that  of  the  actual  crack.  This  is  shown  in  figure  5.4 
which  plots  the  experimental  and  calculated  aspect  ratios  for  test  SB-7. 

As  the  crack  grows,  the  model  predicts  a  slight  decrease  in  closure  at 
the  free  surface,  while  closure  at  the  interior  decre.ases  and  then  increases. 
'Phe  rt-ason  for  the  decrease  in  closure  at  the  free  surface  appears  to  be  the 
increase  in  as  the  crack  length  increases.  'Phe  value  of  K„p  remains 

nearly  constant  as  the  crack  grows.  'Phe  increase  in  interior  clos\ire  at  high 
values  of  a/t  is  probably  due  to  the  inaccuracy  of  the  weight  function 


method  that  was  used  to  calculate  K„p.  Figures  3.5  and  3.6  show  the 
inadequacy  of  the  method  for  largt'  cracks  (a/t  ^  D.it). 

5.2  TESTING  BY  TROH.\: 

I'roha  ,62j  conduct('(l  four  cyclic  hemling  tests  with  surface  Haws  in 
f’MMA  specimens,  'riiese  tests  weri'  |)erforme<i  on  specimens  76.2  mm  wide 
and  19.05  mm  thick,  which  were  ma<‘hined  from  a  ditferent  sheet  of  P.\lM/\ 
than  that  studied  by  Way.  Test  3-11  m.aiiit  ained  a  constant  AK  (d‘ tihO  Kpa 
at  the  int<'rior  (crack  loc.ation  .V)  while  3-11  maintained  the  same  AK 
at  the  free  surface  (cr.ack  loc.alioti  (■).  'I'hese  constant  AK  tests  were 
achieved  by  approjiriatc  slu'dding  of  th<'  applied  cyclic  moment.  .V  con.stant 
cyclic  moment  of  72.8  N-m  was  applie<i  to  3-15,  while  test  3-16  was 
subjected  to  block  loading  (cyclic  moments  of  72.8,  87.3,  72.8,  and  36.1  .\- 
m). 

Troha  observed  phenomena  that  were  not  seen  by  Ray.  Troha  viewed 
what  he  referred  to  as  a  "Type  11"  crack  j63j.  This  type  of  crack,  shown  in 
figure  5.5,  exhibits  closure  at  the  crack  interior,  while  no  contact  occurs  in 
the  central  area  of  the  crack.  Troha  reported  that  this  type  of  crack  had  an 
aspect  ratio  (a/c)  that  ranged  from  0.75  to  0.65  while  a/t  varied  from  0.08‘1 
to  0.156.  As  this  type  of  crack  grew  larger  it  exhibited  the  same  typo  of 
behavior  that  Ray  observed  (no  interior  closure).  This  type  of  crack  was 
referred  to  as  a  "Type  Ill"  crack  by  'I'roha,  which  he  defined  as  a  crack  that 
was  always  open  at  the  interior  (point  A).  In  the  testing  by  Ray,  the  last 
portion  of  the  crack  to  open  was  at  the  free  surface.  For  Troha's  Type  HI 
crack,  the  last  point  to  open  was  about  12  degrees  away  from  the  free 


surface  (figure  5.6). 

f'igtires  0.7  and  5.8  show  oxpc'riinental  and  predicted  closure  for  Troha’s 
constant  load  tc'st  ii-lT).  I’lguri!  5.7  shows  good  agrecuiuMit  for  fref'  surfai'e 
closure.  In  figure  5.8  it  can  be  seen  that  there  was  some  experimentally 
observed  interior  closure  at  small  crack  sizes.  The  model  preilicted  much 
more  interior  closure  than  was  ob.servcd  experimentally.  I'igures  5.9  and 

5.10  show  experimental  and  predicted  closure  for  test  5-16,  which  was 
suhj<a't(.,l  to  blo<-k  loading.  In  figure  5.9  tln're  is  .a  large  increase  m 
experimental  closure  at  a  c/t  ratio  value  of  approximately  0.1.  This  is 
because  the  specimen  experienced  a  fifty  percent  decrease  in  api^lied  moment 
at  this  [)oint.  d'he  predicted  curve  does  not  show  this  dro[)  because 
predicted  the  specimen  would  fail  before  the  load  was  reduced.  Therefore 
another  computer  run  was  made  to  analyze  test  3-16.  In  this  run  crack 
growth  was  not  predicted  by  the  model,  but  instead  forced  to  follow  the 
experimentally  observed  pattern.  The  results  of  this  run  are  shown  in  figure 

5.11  as  the  curve  marked  "Simulated  growth."  Although  this  curve  shows 
the  increase  in  closure  at  c/t=0.4,  it  also  shows  it  decreasing  much  too 
rapidly  as  the  crack  continues  to  grow. 

5.3  TENSION  TESTING: 

Two  cyclic  tension  tests  were  performe<l  with  surface  flawed  plates. 
'ITu;  specimens  were  cut  from  the  same  sheet  of  PMMA  th.at  was  used  by 
Ray,  and  were  oriented  so  that  the  cracks  grew  in  the  same  plane.  It  was 
necessary  to  pre-crack  these  specimens  under  bending  loads  to  avoid  de¬ 
bonding  in  the  grips.  Section  4.4  of  reference  58  explains  how  these 
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specimens  were  prepared. 

As  in  the  bending  tests  conducted  by  Ray,  no  interior  closure  wns 
ob.ser\'ed  for  the  cyclic  tension  loading,  t'onsiderable  closure  ('xistcd  at  the 
fri'e  surface.  'I'his  factor,  along  witii  the  relative  incrtuasc'  in  K  at  the 
interior  caused  the  crack  to  grow  much  more  rapidly  in  tiie  interior  than  at 
the  free  surface.  The  K  at  the  interior  during  bending  was  not  very  high 
since  the  bending  stress  drops  off  quickly  as  one  moves  away  from  the  free 
surface.  Initially,  tlte  free  surfact.'  crack  dimension  c  grew  -eery  slowly  in 
both  specimens  T-2  and  T-3.  .\s  observed  by  Troha  and  Ray,  the  point  of 
maximum  displacement  at  minimum  load  wa.s  internal  rather  than  at  the 
free  surface.  'L'he  last  point  on  the  crack  to  open  was  at  the  free  surface,  as 
Ray  observed,  rather  than  at  12  degrees  inside  the  free  surface  which  Troha 
observed.  Figure  5.12  shows  that  the  model  did  not  correctly  predict  the 
high  free  surface  closure  levels  that  were  observed  experimentally  when  the 
pre-cracking  bending  loads  were  concluded  and  the  tension  loads  initiated. 
In  these  figures  a/t  is  plotted  along  the  x-axis  rather  than  c/t  because  so 
little  growth  occurred  along  the  free  surface.  As  in  the  tests  by  Ray,  the 
model  predicted  significant  interior  closure  when  none  was  observed 
experimentally  (figure  5.12). 

It  is  not  clear  why  so  much  more  closure  is  observed  at  the  free  surface 
than  is  calculated  by  the  model.  The  stress  distribution  under  bending  load 
is  different  from  that  caused  by  tension  loa<rmg  (stri'ss  under  tension  loading 
is  constant  while  stress  under  bending  decreases  as  distanc*'  from  the  free 
surface  increases),  but  it  is  not  apparent  why  this  should  cause  so  much 
closure.  The  high  closure  does  seem  to  be  a  result  of  the  transition  from 
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bending  to  tensile  loading,  since  closure  decreases  as  the  crack  grows  under 
tensile  load. 

5.4  SENSITIVITY  ANALYSIS 

After  the  numerical  model  had  been  completed  and  compared  with  le.sL 
results,  a  series  of  computations  were  made  (see  'Fable  5.2)  to  determine  the 
sensitivity  of  closure  to  several  dilferent  parameters.  The.se  parameters 
were: 

h’low  stress 

*  Ma.ximum  load 

'  R  ratio 

*  Constraint  factor  alpha 

*  Weight  function  reference  case 

Although  the  modulus  of  elasticity  appears  in  some  of  the  equations 
used,  it  has  no  effect  on  the  computed  closure.  Lowering  the  modulus  allows 
greater  displacements  of  the  crack  surfaces  and  the  elements  that  are 
created  are  longer,  but  the  contact  stresses  remain  the  same,  therefore 
closure  remains  the  same.  In  his  literature  survey,  Banarjee  [5]  reports  that 
it  has  been  observed  by  Fuhring  and  Seeger  [63]  that  K^p  shovdd  be 
independent  of  the  elastic  modulus,  but  this  has  not  been  investigated 
experimentally. 

5.4.1  Flow  Stress: 

Yield  stress  is  used  only  to  determine  the  convergence  limit  for  the 
contact  stresses.  When  an  iteration  causes  a  change  in  stress  less  than  two 


percent  of  the  yield  stress  for  all  elements,  the  stresses  are  assumed  to  have 
converged.  In  all  other  calculations  the  flow  stress,  which  is  the  avc.Tage  of 
th<'  yield  and  ultimate  stri'ss,  is  used.  Case  number  1,  shown  in  'ruble  u/J, 
was  bastal  on  Hay's  test  >H-7.  'I'liis  test  had  a  maximum  moment  <d’  70.1 
-Vm,  which  corresponds  to  a  maximum  outer  liber  bending  stress  of  7'J  10 
Kpa  (See  Table  5.1  for  more  information.).  This  experiment  was  modeled 
using  flow  stresses  of  13,790,  27, .580,  41,370,  and  68,950  Kpa  (run  numbers  2, 
3,  1  and  4),  and  the  elfect  of  varying  flow  stress  can  be  seen  in  ligures  5.  HI 
and  5.14.  .\ote  that  the  elfect  of  flow  stress  on  the  closure  calculations 
appears  to  be  negligible.  'I'he  initially  lower  values  of  closure  at  the  free 
surface  for  a  flow  str(>ss  of  13,790  Kpa  are  due  to  an  insuilicitmt  amount  of 
pre-cracking.  When  the  flow  stress  is  lowered,  several  things  happen: 

*  The  plastic  zone  becomes  larger,  and  therefore 
the  effective  crack  length  becomes  larger. 

The  elements  in  the  plastic  zone  and  wake  become  longer. 

*  The  maximum  stress  that  develops  in  the  w’ake  is  smaller, 
but  contact  occurs  over  a  much  longer  distance 

It  is  the  last  elfect  which  causes  the  model  to  initially  predict  lower 
closure  values  for  a  flow  stress  of  13,790  Kpa.  Simulated  pre-cracking  from 
a ^2.82  mm,  c-^3.51  mm  to  a —3.20  mm,  c-^3.58  mm  is  sufficient  for  a  flow 
stress  of  41,370  Kpa,  but  not  enough  for  13,790  Kpa.  'fhis  is  because 
contact  in  the  13,790  Kpa  flow  stress  case  should  extend  behind  the  point 
where  the  .simulated  pre-cracking  began.  Since  the  model  did  not  generate 
any  plastic  wake  in  this  area  (see  figure  5.15),  it  assumes  no  contact  between 
the  crack  surfaces.  This  results  in  a  closure  calculation  that  is  too  low. 
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Banarjee  [5j  reports  that  while  there  is  some  evidence  that  K^p 

incrt'ases  as  yield  stress  decreases,  inatt'rials  with  dilferenl  yield  strengths 
also  h:iV(;  dili'erent  rnicrostrnclural  features.  'Phis  could  result  i.u  a  variation 
in  surface  roughness  whicli  might,  explain  1  Ik-  cliaiigi^  in  c!o'.;iire  associated 
with  yield  stress. 

5.4.2  Maximum  Load: 

According  to  .\ewrnan  19,  chjsurc?  should  <lecr'‘ase  as  the  iiia.xiniuiii 

applied  load  increases.  'Phis  (-ifect  is  most  noticeable  at  higher 

stress  levels  and  lower  R  ratios,  d'aking  c;is(‘  number  1  as  a  b;\selin(',  the 

maxirmim  tnoment  was  increased  to  1.58.2  .N-m  and  237. .5  .\-m.  'Phese  runs 
showf'd  artificially  high  closure  due  to  the  high  crack  growth  r.ate  caused  by 
the  high  loads.  This  caused  the  newly  formed  elements  to  exceed  the  live 
percent  of  the  plastic  zone  criterion.  Runs  that  were  later  made  with 
smaller  load  iticreases  showed  no  significant  change  in  closure.  Banarjee  ;5' 
reports  in  his  literature  survey  that  some  investigators  have  found  closure 
{K„p/K^^^)  to  be  independent  of  vvhile  others  have  observed  a 

decrease  in  closure  as  increases. 

5.4.3  R-Ratio: 

I’igures  .').l()  and  5.17  show  the  elfect  of  R  ratio  on  crack  closure  when 
is  ke[)t  constant  (run  numbers  1,  5,  ti,  7  and  8),  while  figures  5.18  and 
5.19  show  the  R  ratio  effect  when  a  constant  cyclic  moment  is  maintained 
(run  numbers  1,  9,  10,  11  and  12).  Increasing  the  R  ratio  increases  the 
amount  of  closure  because  it  reduces  the  compressive  yielding  in  the  plastic 


wake  at  minimur.i  load.  The  compressive  yielding  shortens  the  elements  in 
the  plastic  wake  which  reduces  the  contact  stresses  and  therefore  closure  as 
well.  It  should  be  noted  that  the  amount  of  pre-cracking  re(i;iir('d  Is 
dependent  upon  the  ii  ratio.  I’or  high  R  ratios,  contact  only  occurs  very 
near  the  crack  tip,  and  very  little  pre-cracking  is  required.  A  computer  run 
with  an  R  ratio  of  -1.  was  made  but  the  results  were  unusable  because  of  .an 
insufficient  amount  of  pre-cracking.  For  Cgures  5.16  and  5.17  the  maximum 
load  for  all  runs  was  kept  constant.  In  figures  5.18  and  5.19  the  cycdic  load 
was  kept  constant  instead  of  the  maximum  load.  These  figures  show  the 
same  trend  that  is  evident  in  figures  5.16  and  5.17.  For  lower  R  values, 
contact  occurs  over  a  larger  area,  and  more  pre-cracking  is  required. 

According  to  an  early  paper  (1971)  by  Fiber  [3],  closure  should  be 
related  to  R  ratio  in  the  following  manner: 

KopA<n,ax=0-5+0.lR+0.dR'  (5.1) 

This  states  that  increasing  the  R  ratio  should  cause  an  increase  in  the 
opening  load.  This  equation  also  indicates  that  K^p/K^^x  should  have  a 
minimum  value  of  0.5  for  positive  R  ratios.  According  to  Banarjee  [5], 
values  of  ranging  from  0.15  to  nearly  1  have  been  reported. 

5.4.4  Constraint  Factor  Alpha: 

The  effect  of  varying  the  constraint  factor  (<t:)  was  also  examined. 
Alpha  is  used  to  simulate  the  degree  of  constraint,  which  affects  the  amount 
of  closure  since  more  closure  occurs  under  plane  stress  conditions  than  under 
pl.ane  strain  conditions.  In  his  analysis  on  aluminum,  Newman  [19]  iis<'d 
(^=1  for  plane  stress  and  ci  -3  for  plane  strain  conditions.  Because  of  its 


higher  poissons  ratio,  o  was  set  equal  to  5  for  plane  strain  conditions  in 
PMMA.  This  is  explained  in  more  detail  in  section  4.1  of  reference  !58j.  Tlie 
degree  of  constraint  around  the  [)erirnet<T  was  d('t,ermined  by  tlie  nnlform 


depth  estimation  method,  ^^hicll  is  ai.-,o  explained  in  section  .4.2  of  I’opc’s 
thesis,  d'est  Sli-7  was  iiuxleled  using  ihrc'e  dilferent  alpha  disLri[)Ulions: 
l.'niform  depth  estimation,  maxirmim  0=^5  (run  nurnbt'r  1) 

*  Uniform  depth  estimation,  maxiimirn  (run  number  14) 

*  •  ■'>  (pl.iiie  rain)  ('Vf'rywiu  re  (run  number  1.4) 

'I'he  results  of  these  c.as'  s  are  displaye^l  in  ligtires  .4.20  and  .4.21  .along  with 
the  experiine-nt  .ally  oi)serv('d  (dosurc?  loads  from  Ray’s  lest  .SB-T.  .\n  attempt 
to  model  <'omplete  pl.ane  stress  around  liie  cr.ack  [)e‘riineter  (  -  -  1 
everywhere)  caused  KOI\  to  exceed  KMAX^  which  resulted  in  a  stoppage  of 
crack  growth  at  A.  This  caused  an  arithmetic  error  to  occur  in  the 
algorithm.  The  program  is  net  able  to  handle  cases  where  the  crack 
completely  stops  growing  at  the  free  surface  or  the  interior. 

Figures  .4.20  and  .5.21  show  that  closure  increases  as  constraint  factor 
alpha  decreases.  This  is  because  lowering  alpha  simulates  conditions  closer 
to  plane  stress  and  allows  greater  yielding.  Allowing  a  maximum  alpha  of  .4 
seems  to  correlate  better  with  the  test  data  than  a  maximum  alpha  of  3. 
Increasing  alpha  increases  the  maximum  stresses  ('ut^)  that  can  occur  in  the 
plastic  zone,  but  it  decreases  the  size  of  the  plastic  zone  (/)).  The  net  result 
is  th.at  le.ss  yielding  occurs,  and  this  results  in  less  (dosiire. 
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6.4.5  Weight  Function  Reference  Case: 
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The  computer  algorithm  utili/es  the  weight  function  method  to 
determine  (.see  .section  :5. 1  of  this  report  for  more  detail).  This  re(iuir(s 

the  u.se  (if  a  reference  cas(;  for  which  the  stress  intensity  factors  and  crack 
displacements  are  known.  'I'o  test  the  ell'ect  of  the  reference  case  on  closures 
calculated  by  the  model,  a  computer  run  vva.s  made  using  a  tension  reference 
case  (case  15)  instead  of  the  bending  reference  case  used  in  all  other  rviiis. 
Comparing  calculations  1  and  1.5  in  (igures  5.22  and  5.23,  it  can  be  seen  that 
the  reference  case  has  little  effect  on  closure.  'I'he  discrepancies  that  exist 
between  the  curves  are  due  to  the  reference  case  approximations  for  the 
surface  crack  displacements  (U^).  The  surface  crack  displacements  are 
necessary  to  evaluate  CUf/Ca  and  in  equations  .3.2  and  3.3.  Initially, 

the  differences  in  closure  loads  between  runs  I  and  15  are  only  two  percent, 
but  this  difference  increases  as  the  crack  becomes  larger.  This  is  consistent 
with  the  results  shown  in  figures  3.5  and  3.6.  These  figures,  which  compare 
the  results  from  the  algorithm  using  the  bending  reference  case  to  the 
Newman-Raju  stress  intensity  factor  solutions,  show  that  the  method  is  less 
accurate  at  larger  crack  sizes.  Mattheck  et  al  [60],  using  a  tension  reference 
case,  also  reported  that  the  method  is  less  accurate  at  higher  a/t  ratios. 

In  summation,  closure  predicted  by  the  computer  model  is  affected  by 
certain  parameters  in  the  following  ways; 

Modulus  of  elasticity  h.as  no  elfect  on  closure 

‘  Flow  stress  has  no  significant  elfect  on  closure 

*  Increasing  the  R  ratio  increa.se3  closure 

*  Increasing  the  constraint  factor  alpha  decreases  closure 

*  Maximum  load  has  no  significant  elfect  on  closure 
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Table  5.1:  Load  levels  and  specimen  dimensions  of  experiments  which  were 
modeled  numerically 
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?• 

I'.xperimcnter 

(Toss 

Section 

Loading 

Load 

Level 

Test 

ID 

R 

K.atio 

Kay 

i()i.i)\-23.n 

milliinetcrs 

bendii\g 

8790  K[)a 
maximum 

SB-1 

0.1 

Ray 

101.6x23.9 

millimeters 

bending 

8790  Kpa 
maximum 

SB-2 

0.1 

Ray 

101.6x23.9 

bending 

10,340  Kpa 

SB-3 

0.1 

millimeters 

maximum 

Ray 

101.6x23  9 

bending 

10.340  Kpa 

SB-6 

0.1 

millimeters 

maximum 

Ray 

101.6x23.9 

Ijending 

7240  Kpa 

SB-7 

0.1 

IS- 

rnilliriieters 

maximum 

Troha 

76.2x19.1 

bending 

AK^=  660 

3-11 

0.035 

millimeters 

Kpa— 

I'roha 

76.2x19.1 

bending 

AK,=  660 

3-14 

0.035 

a 

millimeters 

Kpa— 

Troha 

76.2x19.1 

bending 

16,340  Kpa 

3-15 

0.035 

r-* 

millimeters 

maximum 

T  roha 

76.2x19.1 

bending 

block 

3-16 

0.035 

!■ 

millimeters 

loading 

* 

Pope 

88.9x23.9 

tension 

7000-7930 

T-2 

0.1 

millimeters 

Kpa  maximum 

.  ^. 

Pope 

88.9x23.9 

tension 

7450-6515-5580 

T-3 

0.1 

millimeters 

Kpa  maximum 

Table  5.2:  Summary  of  sensitivity  analysis  runs 
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•1 
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Max  =  5 
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Max  =  5 

41,370 
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Max  =  5 

41,370 

0.3 

Distributed 
Max  =  5 

41,370 

0.5 

Distributed 
Max  =  5 
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— -  5 
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41,370 
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Distributed 

lu-ndiii 

Bendin 


Bending 


Bending 


Bending 


Bending 


Bendin 


Bending 


Ik'ndm 


Bending 


Tension 


{•'igure  5.7:  Free  surface  closure  for  test  .'J-lf)  constant  load 
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F’igure  5.17:  Effect  of  R  ratio  on  closure  at  the  interior 
when  's  kept  constant 
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CHAPTER  B 


SUMMARY 


Th(?  ivo  of  this  (chapter  is  to  sununarizo  th('  ('\[)crin;c:!t  al  aiid 

analytical  results  of  this  project.  Section  B.l  sununarizes  th.e  n-sults 
obtained  from  the  various  th-ough-thickness  Haw  and  part-through  Haw 
experiments  conducted  with  the  transparent  polymer  members.  In  section 
6.2  the  polymer  crack  opening  behavior  is  then  examined  in  context  with 
metal  behavior  reported  by  various  investigators.  The  analytical  results  are 
summarized  in  se'^tion  6.3. 

6.1  SUMMARY  OF  POLYMER  RESULTS: 

Fatigue  crack  closure  experiments  were  conducted  with  through- 
thicknes.s,  corner,  and  surface  flaws  located  in  transparent  polymer 
specimens.  Optical  interference  was  employed  to  determine  the  three- 
dimensional  crack  surface  displacement  field  and  to  characterize  local 
(point-wise)  crack  closure  levels  along  the  crack  perimeter.  In  addition,  bulk 
(average)  crack  opening  loads  obt.ained  for  certain  crack  geometries  v  . 
(.’.MOD  ami  back  face  strain  measurement  techniques,  were  compared  with 
the  point-wise  measure  of  crack  closure  obtained  from  optical  interference. 


The  optical  interference  fringe  patterns  showed  that  a  major  portion  of 
the  crack  surfaces  were  sef)arated  in  tlu?  thick  f)olvrner  sj)ecimeri  at  zero 
loail.  However,  a  rolatlvely  largo  amomil  of  crack  closure  svas  oloir.i  d  at 


the  1  ree  surlace  (plane-stress)  locations.  .\lt  hough  t  he  crack  snrl'acc^  'm  ihe 
specimen  interior  are  prop{)e(l  open  under  zero  load,  th('  l.argi'  hwel  of  free 
surface  closure  alfects  s(jbse<|uent  opening  at  the  interior  location.  For 
example,  in  through-tiiickness  cracked  members,  the  mid-plane  displaciMiient 
prolile  obtained  from  a  numerical  (dastic  .anahsis  agrei'd  well  with  tlie 
e.\p('rimental  results  frorti  o{)tic.al  interference,  provided  that  the  eifect  of 
free  surface  closure  was  accounted  for.  W’lu-n  computing  t  lu;  exj^erimental 
elastic  displacement  profile  at  the  specimen's  mi<l-plane,  it  was  necessary'  to 
select  load  levels  which  caused  a  major  portion  of  the  free  surt'aia;  crack 
surfaces  to  remain  open  (additional  details  are  given  in  reference  ;h0  ). 

Although  crack  tunnelling  was  observed  for  various  flaw  shapes,  a 
quantitative  study  of  the  through-thickness  flaws  showed  that  the 
magnitude  of  crack  tunnelling  was  relatively  invariant  with  crack  length. 
This  behavior  suggests  that  the  crack  growth  rate  along  the  specimen’s  free 
edge  does  not  vary  significantly  from  the  mid-plane  crack  growth  rate, 
although  the  closure  level  at  the  free  surface  is  significantly  higher  than  the 
interior  location.  Thus,  the  effective  cyclic  stress  experienced  at  the  mid¬ 
plane  location  is  not  only  affected  by  the  closure  level  at  that  point,  but  also 
by  the  closure  behavior  along  the  free  surfaces.  This  issue  is  .also 
complicated  by  the  fact  that  the  stress  intensity  factor  for  a  curved  through 
crack  would  be  expected  to  vary  along  the  crack  perimeter. 


Experiments  conducted  with  varying  specimen  thicknesses  showed  that 
crack  closure  was  primarily  confined  to  a  strip  approximately  2-3  rnin  wide 
near  each  free  edge  for  through  Haws  in  I’MMA  test  members.  This 
transition  re  gion  (strip)  was  not  signilica ntly  allected  by  s[)eclmen  thickness 
or  tlse  cyclic  level,  'riius,  a  larger  [a'rcent  age  of  the  specimen’s 

thickness  was  under  this  transition  region  for  the  ’thinner’  members.  Since 
the  difference  in  the  crack  opening  behavior  between  the  free  surface  and 
tfie  rnid-plane  was  less  evidcuit  for  a  thin  member,  thi'  elfective  cyclic  stress 
experienced  at  these  two  locations  did  not  vary  significantly.  As  a  result, 
the  magnitude  of  crack  tunnelling  was  observed  to  decrease  with  a  reduction 
in  specimen  thickness.  (Tack  tunnelling  was  also  observeii  with  corner  Haws 
located  at  holes,  flere  a  comparison  Ix.Tween  the  experimental  crack  front 
with  an  assumed  quarter-elliptic  shape  showed  that  the  actual 
(experimeatal)  crack  front  in  the  specimen  interior  lies  ahead  of  the  elliptic 
shape. 

Comparison  between  the  global  (average)  crack  opening  loads  obtained 
from  CMOD  and  BPS  measurement  techniques  with  the  interferometric 
opening  load  showed  that  for  a  thick  member,  these  bulk  opening  loads  are 
lower  than  the  free  surface  opening  load.  Since  a  thick  member  exhibits  a 
significant  difference  between  the  free  surface  and  the  interior  crack  opening 
behavior,  these  bulk  (global)  crack  opening  loads  are  expected  to  represent 
an  average  quantity  for  the  emtire  specimen  thickness. 
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8.2  CRACK  OPENING  IN  METALS: 

The  objeclive  horo  is  to  compare  aspects  of  the  polymer  results  with  the 
limited  amount  of  iiiformatiou  available  for  experiments  conducted  by 
various  other  investigators  with  metal  members. 

Figures  6. 1-6.3  present  the  crack  opening  loads  obtained  from  various 
measurement  techniques  in  Rene-95  (a  nickel  based  alloy)  specimens.  These 
experiments  were  conducted  with  compact  tension  specimens  by  Dr.  N.  E. 
Ashbaugh,  of  the  University  of  Dayton  Research  Institute  [-19].  Here  the 
metal  specimens  were  subjected  to  constant  cyclic  levels.  Three 

different  techniques,  namely  CMOD,  BPS,  and  interferometric  displacement 
gauge  (IDG)  measurements  were  employed  here  to  give  global  and  local 
crack  opening  load.  The  IDG  method  is  described  in  references  35-38  and 
section  2.2  of  reference  50.  Recall  that  the  CMOD,  and  BFS  readings  reflect 
a  bulk  (average)  crack  opening  load,  while  the  IDG  measurements  provide  a 
local  crack  opening  level  for  a  certain  distance  behind  the  free  surface  crack 


Figures  6. 1-6.3  compare  these  various  measures  of  crack  opening  load. 
The  crack  opening  loads  obtained  from  the  IDG  measurements  were  taken 
at  a  point  7.62  x  10“^  mm  behind  the  crack  tip  on  the  specimen  surface. 
These  figures  compare  the  bulk  crack  opening  (CMOD  and  BFS)  with  local 
free  stirface  (IDG)  crack  opening  for  10,  5,  and  2  mm  thick  specimens. 

Note  from  Figures  6.1  and  6.2  that  for  relatively  thick  members,  the 
bulk  opening  loads  fall  significantly  below  the  free  surface  opening  loads. 
This  difference  in  crack  opening  was  also  observed  with  the  polymer 
specimens,  and  suggests  that  these  thick  metal  members  also  exhibit 


significant  three-dimensional  variation  in  crack  opening.  On  the  other  hand, 
Figure  ti.H  shows  that  for  a  thin  member  (2  mm  thick),  the  difTerence 
between  the  bulk  and  free  surface  opening  is  not  significant,  suggc'stiiig  .an 
uniform  crack  o[)('ning  behavior  in  the  thin  members. 

I'his  variation  in  crack  opening  between  ’thick’  and  ’thin’  members  was 
also  observed  with  the  PMMA  through-flaw  experiments,  where  the  bulk 
(average)  measure  of  crack  opening  load  for  a  thick  member  was  found  to  be 
significantly  less  than  the  interferometric  free  surface  opening  load. 
Furthermore,  interferometric  results  obtained  from  polymer  experiments 
conducted  with  varying  specimen  thickness  showed  that  crack  opening  for  a 
thin  (.a  mm  thick)  FM\L\.  specimen  was  uniform  across  the  crack  plane, 
while  a  25  mm  thick  inember  indicated  a  large  displacement  gradient 
between  the  mid-piane  (interior)  and  the  free  edges.  Thus  the  three- 
dimensional  variation  in  crack  opening  decreases  with  a  reduction  in 
specimen  thickness. 

Since  the  difference  in  crack  opening  between  the  free  surface  and  the 
mid-plane  was  le.ss  evident  for  a  thin  member,  the  thinner  specimens 
exhibited  significantly  less  crack  tunnelling  than  the  thick  members.  This 
decrease  in  crack  tunnelling  with  specimen  thickness  was  also  observed  in 
HT80  steel  specimens  subjected  to  overloads  [9]. 

(Jrack  closure  experiments  conducted  with  otner  metal  members  in 
references  j65-ti7|  employing  electron  fractography  showed  a  dilfcrence  in  the 
crack  opening  behavior  between  the  free  surface  and  the  mid-plane 
(interior).  Here  the  specimen's  fracture  surfaces  were  marked  with  a  special 
type  of  !oad-se<pience,  and  fractographic  studies  clearly  indicated  that  the 


free  surface  crack  opening  load  along  the  free  surface  of  a  thick  member  was 
signiiicantly  larger  than  the  mid-plane  crack  opening  value.  Other 
ext)eriiiients  conducted  with  a  push-rod  displacement  gage  in  metal 
sfx’ciimuis  ;U)-hlj  also  show(nl  similar  crack  opening  behavior,  wheri;  the 
crack  first  opens  in  the  interior  followed  by  subsequent  opening  along  the 
free  edges. 

Recently  a  three-dimensional  elastic  plastic  cyclic  finite  element  analysis 
of  a  straight-through  crack  front  |l8j  predicted  the  variation  ifi  fatigue 
crack  opening  load  along  the  flaw  perimeter.  Here  a  straight-through  crack 
located  in  a  finite  thickness  plate  was  simulated  with  a  tensile  cyclic  load 
range  between  2.5  and  25%  of  the  material’s  yield  stress.  The  elastic-plastic 
finite  element  analysis  was  carried  out  for  10  cycles  of  load.  The  crack  was 
allowed  to  extend  one  element  size  (0.03  mm)  each  cycle,  and  the  opening 
loads  were  computed  after  each  load  cycle.  In  this  model,  as  before,  the 
crack  first  opened  in  the  plane  strain  interior  followed  by  crack  tip  opening 
at  the  free  surface  locations  as  schematically  shown  in  Figure  6.4.  Note  here 
that  under  an  applied  load  of  0.1  the  interior  portion  of  the  crack 

starts  to  open.  As  the  applied  load  is  increased  to  0.36  a  major 

portion  of  the  original  interior  crack  tip  remains  open,  while  the  the  free 
surface  crack  faces  are  still  closed.  A  load  level  greater  than  0,56  is 

required  to  separate  the  crack  faces  completely  along  the  free  edges.  Recall 
that  the  through-thickness  flaw  experiments  also  showed  that  the  crack  first 
opened  in  the  specimen  interior  followed  by  free  surface  opening  at  a 
significantly  larger  load. 


Note  in  Figure  6.4  that  unlike  the  polymer  results,  where  a  major 
portion  of  the  interior  crack  perimeter  remains  open  under  zero  load,  the 
numerica!  analysis  in  r<-‘fer(!nce  shows  the  <Tack  I'aiu's  in  the  specimen 
interior  to  he  in  contact  at  zero  lo.a<l.  However,  the  contact  stress(!s  at  tiie 
specimen’s  mid-plane  were  fouml  to  l)e  signilicantly  less  than  those 
computed  at  the  plane  stress  free  edges.  Furthermore,  the  numerical  study 
does  not  include  the  roughness  induced  crack  closure  phenomenon,  where  the 
crack  surface  roughness  keeps  the-  crack  propixnl  opem  under  zt^ro  lo:id. 
Thus,  it  is  believed  here  that  the  initial  residual  crack  opening  in  the 
polymer  specimen  interior  may  be  due  to  a  combination  of  crack  surface 
roughness  and  low  ’contact  stress’  level  observed  at  the  plane  strain  interior. 
.\lthough  similar  crack  surface  roughness  exists  along  the  free  surfaces,  the 
’large’  contact  stresses  force  the  free  edges  to  come  in  close  contact,  and 
prevent  the  crack  from  opening  under  a  relatively  low  load  at  these  plane 
stress  locations. 

Although  a  complete  three-dimensional  experimental  investigation  of 
crack  opening  in  metal  members  may  not  be  po.ssible,  the  indirect  evidence 
presented  above,  such  as  crack  tunnelling  and  fracture  surface  appearance, 
clearly  suggest  a  three-dimensional  variation  in  crack  opening  along  the  Haw 
perimeter  of  a  metal  specimen.  Although  it  is  not  known  whether  the 
polymer  results  obtained  here  describe  the  crack  opening  in  metal  members 
quantitatively,  it  is  believed  that  the  crack  opening  behavior  in  polymer 
specitneris  may  rellect  the  corresponding  behavior  in  the  metal  members 
qualitatively. 
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6.3  SUMMARY  OF  ANALYTICAL  RESULTS: 


Comparison  of  the  closure  calculations  performed  here  with  the  hencling 
tests  by  Hay  and  the  tension  tests  performed  by  Hope  reveals  a  serious 
deliciency  in  the  model,  d’lu'  model  [)redicts  contact  betwu'en  the  crack 
surfaces  at  the  specimen  interior  where  no  contact  is  observed  (see  figure 
6.5).  In  Troha’s  experiments,  some  interior  contact  was  observed,  but  only 
for  small  a/t  ratios  (less  than  0.156).  His  experiments  were  performed  with 
a  thinner  sheet  of  P\1.\I.\  (0.75  versus  0.9-1  inches)  and  at  a  slightly  lowa-r  H 
ratio  (0.035  versus  O.l) 

Flock’s  paper  (20j  on  plasticity-induced  plane  strain  closure  may  explain 
this  absence  of  interior  contact  observed  in  the  experiments.  In  liis  study,  a 
center  cracked  panel  and  an  edge  cracked  bend  specimen  were  examined 
with  a  two-dimensional  elastic-perfectly  plastic  finite  element  code.  Larsson 
and  Carlsson  [68j  showed  the  importance  of  the  T-stress  (the  non-singular 
stress  parallel  to  the  crack  plane)  for  stationary  cracks,  and  Fleck  examined 
its  effect  on  a  growing  crack  under  plane  strain  conditions.  Along  slice  0-6 
in  figure  3.3,  the  surface  crack  may  be  assumed  to  act  like  an  edge  crack. 
The  T-stress  for  an  edge  crack  in  bending  is  tensile,  which  results  in  a 
plastic  zone  at  maximum  load  like  that  in  figure  6.6.  Here  a^,  is  the  crack 
length  at  which  the  finite  element  analysis  was  initiated,  and  af  is  the  final 
crack  size.  The  plastic  zones  shown  in  figure  6.6  occur  when  the  crack 
length  is  n.f.  The  portion  of  the  plastic  zone  along  the  crack  flank  yields  in 
tension  in  the  x  direction,  d'his  counteracts  the  earlier  yielding  that 
occurred  in  the  plastic  zone  ahead  of  the  crack  tip  and  reduces  the  size  of 
the  plastic  wake  elements.  The  net  result  is  less  closure  since  there  is  less 


material  to  wedge  the  crack  open. 

All  experiments  reviewed  in  this  report  showed  that  the  maximum 
displacenuuit  did  not  occur  at  the  crack  mouth,  as  might  be  ex[)e(-t<,'d.  but 
occurred  at  the  crack  interior,  h'h'ck’s  pap<'r  also  sheds  some  light  on  this 
phenomenon.  A  crack  under  plane  stress  closes  continuously,  which  means 
that  contact  first  occurs  immediately  behind  the  crack  tip  and  then  proceeds 
back  toward  the  crack  mouth.  In  a  plane  strain  situation,  discontinuous 
closure  occtirs  after  the  crack  has  grown  a  short  distance  from  its  initial  Haw 
size.  Discontinuous  closure  means  that  contact  first  occurs  immediatcdy 
behind  the  crack  tip  and  next  occurs  far  behind  the  crack  tip  as  shown  in 
figure  6.7.  In  this  figure,  the  gap  between  the  crack  surfaces  at  minimum 
load,  ('’rnin,  divided  by  the  quantity  is  plotted  along  the  ordinate. 

The  distance  behind  the  crack  tip,  x,  divided  by  plotted  along 

the  abscissa.  Figure  6.7  shows  the  displacement  profiles  for  both  center 
cracked  panel  and  bend  specimens.  Notice  that  the  point  of  maximum 
displacement  is  much  closer  to  the  crack  tip  than  the  crack  mouth.  This  is 
the  same  type  of  behavior  that  has  been  observed  in  the  specimen  interior 
for  the  surface  crack  experiments. 

Ray  conducted  experiments  at  three  different  stress  levels.  As  expected, 
the  medium  stress  level  exhibited  higher  free  surface  opening  loads  than  the 
lower  stress  level  tests.  The  highest  stressed  tests  had  lower  free  surface 
opening  loads  than  the  medium  stressed  tests.  The  computer  model  does  not 
predict  this  trend.  As  the  applied  load  is  increased,  the  model  will  always 
predict  higher  opening  loads.  This  discrepancy  may  be  due  to  the  effect  of 
the  T-stress.  Fleck  reported  that  for  plane  strain  conditions,  no  plasticity- 


induced  closure  would  occur  if  the  ratio  greater  than  a  critical 

value  between  0.035  and  0.070. 

In  Troha’s  tests,  he  notic(!d  that  as  load  w;is  applied,  the  last  point  <;n 
the  crack  perimeter  to  open  was  tiot  at  the  free  surface  but  about  twelve 
degrees  inside  the  free  surface.  Elastic  plastic  Unite  element  analysis  by 
Trantina  [69j  on  semicircular  surface  cracks  showed  that  the  highest  stress 
occurred  at  about  fifteen  degrees  from  the  free  surface,  d’his  was  due  to  the 
redistribution  in  load,  and  may  be  responsible  for  the  opening  bidiavior  that 
Troha  observed.  His  specimens  were  thinner  than  those  of  the  other 
experimenters,  which  may  be  why  this  phenomena  was  not  observed  by 
others. 

The  computational  model  used  in  this  project  is  based  on  the 
plasticity-induced  closure  mechanism  only.  The  roughness  induced 
mechanism,  which  is  usually  significant  only  for  small  crack  or  near 
threshold  cases,  is  not  accounted  for  in  this  model.  In  most  cases  analyzed 
in  this  paper,  the  predicted  closure  is  greater  than  the  actual  closure,  so  the 
inclusion  of  a  roughness  induced  mechanism  would  not  reduce  these 
discrepancies.  Visual  inspection  of  the  fracture  surfaces  of  specimen  T-2 
showed  smooth  transparent  surfaces.  Also  the  interference  fringes  appeared 
to  be  reasonably  distinct.  In  the  case  of  polycarbonate  [7],  the  interference 
fringes  cannot  be  seen  when  the  crack  is  grown  at  high  values  of  AK,  which 
causes  the  crack  surfaces  to  be  become  rough. 

From  this  project  the  following  conclusions  can  be  drawn: 

1.  The  free  surface  closure  values  predicted  by  the  model 
correlate  well  with  most  of  the  experiments  examined. 


2.  The  model  predicts  closure  in  the  interior  where  none 
or  very  little  is  observed. 

i}.  At  ininiimiin  load  the  inaxiinuin  disj)laceinent  is  ob.s<'rved 
to  occur  at  an  internal  location,  whilt-  the  model 
assumes  it  occurs  at  the  crack  mouth. 

4.  Experiments  show  that  K^p  may  decrease  when  load 
is  increased.  The  model  predicts  K„p  always  increases 
as  the  load  increases. 

The  basic  method  of  the  model  appears  to  be  valid,  but  the  equations 
used  to  calculate  yielding  and  displacement  seem  to  be  inadequate.  .Since 
the  observed  contact  is  very  different  from  that  predicted  by  the  model 
(figure  6.5),  either  the  calculated  displacements  of  the  crack  surface  or  the 
lengths  of  the  plastic  wake  elements  are  incorrect.  In  the  current  model,  the 
effect  of  yielding  along  the  crack  faces  in  the  plastic  wake  is  ignored.  This 
could  cause  elements  in  the  plastic  wake  to  be  too  large,  which  would  result 
in  erroneously  high  values  of  Kyp.  In  determining  the  displacement  of  the 
surface  crack  under  remote  loading,  it  is  assumed  that  the  crack  displaces 
along  various  slices  (figure  3.1)  like  a  center  cracked  panel.  This  results  in 
the  maximum  displacement  occurring  at  the  crack  origin.  Under  minimum 
load  it  is  observed  that  the  maximum  displacement  occurs  well  inside  the 
crack.  This  indicates  that  the  displacement  of  the  surface  crack  is  much 
more  complex  than  the  model  assumes. 

Several  recommendations  can  be  made  for  future  research  work  on 
surface  crack  closure  modeling; 

1.  A  three-dimensional  elastic-plastic  finite  element  analysis  of  a 
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growing  surface  crack  should  be  performed.  This  type  of  analysis 
has  already  been  done  for  an  edge  crack  by  Chermahini  [I8j,  and  if 
pc-rforiiuai  for  a  surface  crack,  would  give  a  better  understanding 
of  tlie  contact  stresses  oceurring  in  the  plastic  wake. 

2.  A  beltt'r  elastic  displacement  efjuation  for  the  surface  crack 
needs  to  be  developed.  The  method  used  in  the  present  model 
incorrectly  predicts  the  location  of  the  maximum  displacement. 

More  tesfing  should  be  pi  rformed  to  det(;rmine  the  effect  of 
specimen  thickness  on  stirftice  crack  beh.avior.  Some  differences 
in  behavior  have  imen  noted  between  Ray’s  0.94  and  Troha’s  0.7.5 
inch  thick  specimens. 

4.  More  testing  should  be  performed  to  determine  the  effect  of  load 
level  for  constant  amplitude  loading.  Testing  by  Ray  has  indicated 
that  Kqp  can  decrease  or  increase  as  applied  load  is  increased. 
Accomplishing  these  four  items  should  lead  to  a  better  understanding  of 
closure  in  surface  cracks  and  facilitate  improved  closure  modeling. 
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Figure  6.6:  Plastic  zone  for  an  edge  crack  {reference  [20]) 


•"0  noae  to  Cose 

and  penultimate 
ncde  to  open,  for  ^ 
3end  ,CCP' 


-C3 


initial  pcsit  on 
of  CrGCi<  tip 
(0:0.) 


Crccl<  tip 


Figure  6.7:  Discontinous  closure  (reference  [20]) 
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